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ABSTRACT 

A proposed ship propulsion system which incorporates 
superconducting electric machines as the transmission 
system between the prime mover and the propeller is 
described. The propulsion system employs gas turbine 
prime movers, synchronous generators and synchronous 
motors with superconducting field windings, switch gear 
with a cycloconverter, variable frequency, power controller 
between the generators and the motors. The proposed 
system in the DD963 destroyer, which has a gas turbine 
propulsion system driving controllable pitch propellers 
through reduction gears. The resulting ship is compared 
with the original on the basis of weight and volume. A 
smaller ship with an identical payload but a smaller pro- 
pulsion system is constructed to take advantage of the 
weight and volume savings which are a result of using 
superconducting electric machinery. The smaller ship is 
compared with the original DD963 on the basis of weight, 
volume, effeciency and cost ceiling for the superconduct- 
ing electric propulsion system. 

The proposed superconducting motors and generators 
are modeled mathematically and simulated on a digital 
computer. Components for the motors and generators are 
designed to determine their individual characteristics 
and their interactions with other elements of the machines. 
The design analysis of the superconducting machines 
indicates they will be very small and lightweight. 

Final comparison of the proposed and existing ships 
shows a 14 # reduction in ship displacement, a 9 % reduc- 
tion in total volume, a 1 7# reduction in fuel carried and 
a propulsion system that is 50# lighter and requires 33 ^ 
less volume. 

Thesis Supervisor: J.L. Kirtley, Jr. 

Title: Associate Professor of Electrical Engineering 

Thesis Reader: J.W. Devanney 

Title: Associate Professor of Marine Systems 
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NAVAL ARCHITECTURE DEFINITIONS 



SYMBOL UNITS 



DEFINITION 



I. Gross Characteristics 



A 


(tons) 


Fuel Load Displacement 


V 


(ft 3 ) 


Total internal volume 


L or LBP 


(ft) 


Length of hull at waterline 


B 


(ft) 


Beam(width) of hull 


d avg 


(ft) 


Average depth or average 
height of main deck 


T 


(ft) 


Draft of hull 


VCG 


(ft) 


Vertical Center of Gravity 
of ship 


C o 




Prismatic coefficient. The 



percentage of a prism, which 
is the same L,B,T as the 
underwater hull , actually 
filled by the underwater 
hull. For a Box of L,B,and 
T units, C = 1. 



C x 




Midship section coefficient. 

The percentage of a rectangle, 
which is the same B and T as 
the underwater hull, actually 
filled by the midship hull 
cross section. For a rectangle 

C x = 1 - 


GM/B 




Measure of ship stability 
(Resistance to rolling) 


V 

s 


(knots) 


Maximum sustained speed 


R$V„ 

e 


(Naut.Mi.@Kts) 


Range at cruising speed 


SHP 


(hp) 


Main propulsion horsepower 
rating 


SFC 


( lbs/hp-hr) 


Specific fuel consumption in 
pounds per horsepower-hour 


M 


(M) 


Total manning complement 
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SYMBOL 



UNITS 



DEFINITION 



II. Weight Fractions* 

WTGPl/ a 

W MB/A 
WTGP2/ a 

W OPS/a 

W PERS/a 
W PAY/ a 



Structural weight fraction 

Machinery Box weight fraction 

Propulsion system weight 
fraction 

Ships operations weight 
fraction 

Personnel weight fraction 
Payload weight fraction 



III. Volume Fractions 

VOL MB/v 
VOL OPS/v 

VOL PERS/v 
VOL PAY/v 

* weight and volume fractions 

IV. Specific Ratios 



Machinery Box volume fraction 

Ships operations volume 
fraction 

Personnel volume fraction 
Payload volume fraction 
are non-dimensional 



WTGP2/SHP (lbs/SHP) Propulsion system specific 

weight 

VOL MB/SHP ( ft 3 /SHP) Machinery Box specific 

volume 

VOL HAB/M (ft 3 /MAN) Personnel specific volume 

V. Densities 

DISPLACEMENTAOLUME (lbs/ft 3 ) Total ship density 

VI. . BSC I Weight Groups * 

BSCI weight groups are a breakdown of ship weights by ship 
systems as listed below: 

WTGPl (tons) Hull Structure 
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SYMBOL 


UNITS 


DEFINITION 


WTGP2 


(tons) 


Propulsion System 


WTGP3 


(tons) 


Electric Plant 


WTGP4 


(tons) 


Communication and Control 


WTGP5 


(tons) 


Auxiliary Systems 


WTGP6 


(tons) 


Outfit and Furnishings 


WTGP7 


(tons) 


Armament 


WTGP8 


(tons) 


Loads 



* Complete listing of the contents of each weight group 
are contained in Appendix I. 
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INTRODUCTION 



Superconducting electrical machines can now be 
considered "state-of-the-art" for marine applications. 
This paper investigates the impact of converting a DD963 
gas turbine driven propulsion plant from mechanical power 
transmission to superconducting electrical power trans- 
mission. The DD963 is a 7885 ton, twin screw, 30+- knot 
destroyer with a 6000 mile endurance range. The propul- 
sion plant prime movers are four 20,000 horsepower gas 
turbines. The mechanical drive propulsion system of the 
DD963 requires two gas turbines driving the reduction 
gear for each shaft. The physical connection of a gas 
turbine to a reduction gear limits a gas turbine to 
driving that one shaft only. In a electric drive pro- 
pulsion system any one gas turbine can drive either or 
both shafts at the same time. At the cruising speed of 
20 knots, the mechanical system required two gas turbines 
(one for each shaft) to be in operation. At this low 
power level, each turbine is operating in a very unecono- 
mical off-optimum fuel consumption performance mode. A 
single gas turbine providing the power for a 20 knot 
cruising speed and using electric drive to transmit 
this power to both shafts, operates at a more economical 
performance power level. For this reason, an electric 
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drive requires much less fuel than the mechanical drive 
system. The criteria for comparing the electric drive 
ship to the original mechanical drive DD963 is the ships 
must perform the same military mission without changing 
the speed and endurance characteristics. 

Figure 1 contains a major weight breakdown for the 
DD963 baseline ship, a four gas turbine electric drive 
ship, and a three gas turbine electric drive ship. When 
a straight conversion of the four gas turbine driven 
DD963 to electric drive is made a propulsion plant 
weight and fuel weight savings of 485.1 tons is realized. 
Since the DD963 is a weight limited ship (no weight margin 
to spare) , this is a dramatic improvement in the overall 
ship characteristic. 

The DD963 was originally designed with a large 
excess margin of volume, which is now even greater with 
the 485 ton reduction in displacement. The next logical 
step was to take full advantage of the propulsion plant 
and fuel weight reductions and reduce the overall size 
of the ship. This further reduced the fuel weight* less 
fuel required to drive a smaller lighter weight ship. 

A smaller and lighter ship also requires less installed 
horsepower, which permitted the reduction of the propul- 
sion plant to three 20,000 horsepower gas turbines. An 
overall savings due to hull, propulsion plant and fuel 
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weight reductions in a three gas turbined powered ship 
is 1084 tons for a 14# weight reduction over the base- 
line ship. 

The three- engined ship has the same mission perform- 
ance capabilities as the larger four-engined ship. The 
only differences between the two ships is the three- 
engined ship is cheaper to build, maintain and operate. 



DD 963 4 ENGINED 3 ENGINED 

BASELINE ELECTRIC DRIVE ELECTRIC DRIVE 



Hull Structure 


3137.1 


3105.6 


2757.7 


Propulsion 


789.2 


504.1 


401.2 


Electric Plant 


296.8 


296.8 


275.9 


Command & Control 


250.3 


250.3 


250.3 


Auxiliary Systems 


739.8 


739.8 


735.5 


Outfit & Furnishings 


454.3 


454.3 


445.9 


Armament 


159.2 


159.2 


159.2 


Margin 


100 


100 


100 


Fuel 


1606 


1404 


133 ^ 


Loads 


353 


353 


341 


Full Load 


Displacement 


7885 


7366.4 


6800.7 


% Weight Saving from 


Baseline 


0 


6# 


14# 


MAJOR 


WEIGHT 


BREAKDOWN 





FOR THE DD 963 BASELINE SHIP 
AND TWO ALTERNATIVE ELECTRIC DRIVE SHIPS 



Figure 1 
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In comparison to other propulsion systems, the super- 
conduction electric machinery offers a truly compatable 
transmission system with which to take advantage of the 
high-speed, compact and highly maintainable gas turbine, 
without the use of gears or controllable reversible pitch 
propellers for reversing. The arrangement flexibility 
inherent to electrical propulsion systems can now be 
realized without the weight penalty associated with con- 
ventional electrical machines. The cost of installing 
superconducting machines onboard a ship appears to be 
feasible when compared to the potential benefits associated 
with this use. 
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CHAPTER 1 



ELECTRICAL MACHINERY FOR SHIP PROPULSION 



1.1 Background 

Interest in electrical propulsion systems has been 
generated by their inherent ability to provide speed 
reduction between a high-speed, efficient, lightweight 
prime mover and a much slower, efficient propeller. The 
primary advantages of such an installation is the flex- 
ibility of design and arrangement of the machinery plant 
and the flexibility of control. 

While electrical transmission systems have a number 
of advantages, the primary reasons for a lack of wide- 
spread use have been: 

Higher acquisition cost than mechanical drive 
alternatives. 

Greater weight and volume requirements than mechani- 
cal drives. 

Higher transmission losses overall, resulting in 
a lower total system efficiency and a higher fuel 
usage than mechanical drives. 

With the development of superconducting electrical 
machinery for shipboard use, a great savings in weight, 
cost and volume may now be obtainable. Superconducting 
machinery provides all of the advantages of electrical 
propulsion with the disadvantages of large weight and 
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volume being greatly reduced. The reductions in size and 
volume can have a dramatic effect upon the ship design 
by eliminating the need for the large machinery box 
required by the current mechanical drive propulsion 
systems. The result will be an increase in available 
volume in a highly desirable location, which can then 
be used for other shipboard functions. Another alternative 
result of a volume and weight saving propulsion system 
may be a smaller and less expensive ship that performs 
the same missions as the larger volume ship with a 
mechanical drive system. 

In the case where the power plant requires the use 
of several prime movers, the electric drive provides 
an efficient method of coupling these units to the 
propeller without the use of mechanical clutches or 
couplings. The electric drive system can be arranged 
in such a manner that the ship operating at less than 
full power will require only a minimum number of prime 
movers to be in service. (Chapter 2 contains the propulsion 
plant operating characteristics.) This contributes to 
greater fuel economy and provides down time for routine 
maintenance on idle propulsion units. 

Since some prime movers, such as gas turbines, are 
unidirectional machines; an electric drive can produce 
the required reverse rotation of the propeller by relatively 
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simple controls. This eliminates the requirement of the 
unidirectional prime movers for controllable and reversable- 
pitch propellers (CRP) to provide the desired reverse 
rotation. Switching to fixed-pitch propellers eliminates 
several disadvantages of the CRP: the extensive hydraulic 
control system, the maximum upper loading limit of a single 
propeller to 40,000 shaft horsepower, and the reduced 
effeciency of a CRP propeller over a fixed pitch propeller 
of the same given size and characteristics. 

Electric propulsion systems are classified as 
either direct-current (dc) or alternating- current (ac). 
Electric systems can further be defined by the type of 
prime mover involved, such as diesel engines or gas 
turbines. Traditionally dc drives have been desirable 
because they provide more rapid and continuous control 
of the propeller speed for excellent maneuverablility. 

As a result, superconductors were first applied to dc 
machines, thus increasing the maximum practical power 
level while allowing for small, light weight machinery. 
Superconducting dc machines require that full electric 
power must be carried onto the rotor at high current 
and low voltage. This presents significant current 
collection problems in high-power machines and constitutes 
a distinct disadvantage . 2) ( 3) 

The high-power electrical connections for an ac 
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machine can be made directly to a stationary armature 
winding. They are of high voltage and low current when 
compared to a dc machine and are not subject to the high 
current collection problems of the dc machine. However, 
a high-power ac machine will require a rotating super- 
conducting field winding of low current and voltage. 
High-power ac synchronous machines have been proven 
practical by experimentation at MIT^ 2 ^^) (*0 an( j elsewhere. 
A great deal of work has also gone into proving that ac 
machines are suited for shipboard use.^^^^ 

The superconducting synchronous motor and generator 
under investigation in this paper belongs to the class 
of machines referred to as cold shield superconducting 
machines. (see Fig. 1.2) The principal functions of 
the cold shield (referred to as cryogenic shield or shield) 
cure to shield the superconducting field winding from 
alternating magnetic fields and to prevent heat transfer 
in the form of thermal radiation. The rotor itself 

is held at or near 4.2°K while the shield would operate 
at a temperature of about 20°K. The damper shield (called 
the damper) operates at approximately room temperature 
and serves as an electro-mechanical damper and as a shield 
for time-varying fields. In the event of a fault, the 
damper absorbs strong crushing and torque loads. These 
loads are so strong that the damper must be strong, 
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CROSS SECTION OF SUPERCONDUCTING MACHINE 



s 

Figure 1.1 
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Figure 1.2 

SU PERCONDUCTING MACHINE ELEMENTS 



consequently thick, relative to the shield. ^ 

The rotor cross section consists of eight elements. 

(see Fig. 1.1) Arranged from inside to outside, they 
ares torque tube, superconducting field winding, field 
hoop support, shield insulation gap, shield, shield 
support, damper insulation gap, and damper. Between the 
armature and the damper is a gap which holds a vacuum to 
reduce the rotor windage loss and to provide thermal 
insulation. 

The outer shell outside the armature (called the stator 
core) provides a uniform boundary condition and confines 
the magnetic field to the machine. This shell is of 
laminated iron to reduce eddy current losses. The iron 
in this shield is soft and must be surrounded by an outer 
shell to protect it. This steel outer shell acts as a 
frame and serves as a structural support for the entire 
machine. Fig. 1.1 

The torque tube is actually a cylinder that serves 
as a cryogenic thermal distance piece and supports the 
field winding. The torque tube must be thick enough to 
withstand the torque that is imparted to the cylinder by 
the magnetic flux field. The torque tube is sized for 
normal torque, based on the machine's full power rating. 
Stainless steel is used for the construction of the torque 
tube. 
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A superconducting field winding is the heart of the 
superconducting machine. This superconducting field 
produces an intense magnetic field without the use oft 
heavy and bulky ferrous material, electrical dissipation, 
and negligible electric power losses. When the temperature 
is reduced below a critical value (approximately 5°K) the 
superconductors support very high currents without 
resistance losses. Current in the stationary armature 
winding interacts with the large flux wave generated by 
the field winding and pulls the rotor around at synchro- 
nous speed in the same manner as a conventional machine. 

A conventional electric motor or generator operating at 
room temperature requires a heavy iron core to produce 
the magnetic fields necessary for proper operation. At 
these temperatures, there is electrical dissipation and 
power losses. To overcome these losses, the machine must 
be made even larger. The size and weight of the iron 
core then controls the size of the conventional machine. 
Conventional machines become big and heavy when compared 
to superconducting machines of the same size. It is in 
this manner that a much smaller superconducting machine 
develops the same horsepower as a much larger conventional 

machine. <3H 5) (7) (9) 

1.2 Baseline Ship 

A particular propulsion system cannot be judged 
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"as good as" dr "better than" another propulsion system 
unless "the systems are compared for a particular mission. 

For a Navy ship the primary mission of the propulsion 
system is to move a given pay load over a given distance 
at a required speed. Simple comparison of one propulsion 
system to another is not sufficient. The overall impact 
of each system on a given ship must be determined. 

Since, in most cases, the propulsion system has the largest 
single impact on the total ship, the propulsion plant 
becomes a major factor in final ship size, cost and pay 
load.^®^ The sizeable volume and weight required for 
a propulsion plant and its fuel over shadows all other 
volume and weight requirements in a Navy ship. Not only 
is the magnitude large, but a feature unique to the 
propulsion system requires space which cannot be split 
up or scattered throughtout the ship, and in most cases 
it occupies the prime space in the ship.^^' For the 
above reasons, a ship with a conventional propulsion 
system was needed for propulsion system comparison* Such 
a baseline ship was found in the DD963 class ship. The 
selection was based on the following primary consideration: 
The best lightweight prime movers for the superconducting 
system are gas turbines {therefore, the baseline ship 
should be gas turbine powered. 

The selection of the DD963 as "the baseline ship was 
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governed by the fact that it is the only totally gas 

turbine powered ship in the U.S. Navy about which a good 

deal of specifications and information has been published. 

The basic characteristics of the DD 963 are contained in 

Table 1.1. ( 12 )d3) Th e p 0wer plant weights by BSCI 

( 12 } 

subgroups are shown in Table 1.2. ' The discrepancy 

between weight group 2 in Table 1.1 and the propulsions 
plant total weight in Table 1.2 can probably be attribu- 
ted to two different authors assigning individual aux- 
iliary equipment weights to different weight groups. 

For a comparison of propulsion systems to be 
based on a computer synthesis model, the volume and 
weight associated with the remaining ship functions 
must be held constant. The items to be held constant 
are BSCI weight groups 3 *^>5 .6 and 7* (see Table 1.1 for 
definition of BSCI weight groups) . The full load loads 
will be allowed to change only to the extent required 
by changes in fuel dictated by the respective power 
plants. Weight group 1 will change by the amount 
required to compensate for different propulsion plant 
sizes and the different fuel requirements associated with 
each propulsion plant. With the above criteria being 
observed throughout the computer simulations, the re- 
sultant weight and volume changes will be due only to 
propulsion plant changes. The selection of the "best" 
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Length 

Beam 

Draft 

Speed 

Displacement 

Crew 

Propulsion 



Electrical 

Power 

Armament 



Sensors 



BSCI GROUP 
1 
2 

3 

4 

5 

6 

7 



529 ft. water line, 563.3 ft. overall 
55 ft. 

29 ft. (navigational) 19 ft. (hull) 

30 + knots 

Approximately 78OO Long Tons (fully loaded) 
Approximately 18 officers 232 enlisted men 
4 LM2500 Gas turbine engines 
80,000 shp 

2 shafts, 2 Controllable Pitch Propellers 



3 Gas turbine driven generators 
2-5” * 54 caliber guns 

1 ASROC 8 Tube Launcher 

2 ASW Torpedo Mounts-Triple barrel 
Fire Control, Surface Search and Air 
Search Radars, Long Range Sonar 



DESCRIPTION 
Hull Structure 
Propulsion 
Electric Plant 
Communications & Control 
Auxiliary Systems 
Outfit & Furnishings 
Armament 

Light Ship (W/0 Margin) 
Margin 

Light Ship (with Margin) 

Full Load Loads 

Full Load Displacement 



WT (LONG' TONS) 
3105 . 5 ^ 
759.93 
293.81 
354.17 
722.98 
452.01 
152.16 
5829.61 
89.87 
5919.48 
I865 . 66 
7785.14 



Table 1.1 
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BSCI NO. 


SUBGROUP DESCRIPTION 


WT(L0NG TONS) 


201 


Propulsion Units 


244.14 


203 


Shafting, Bearings & Propellers 


253.15 


204 


Combustion Air Supply Systems 


58.34 


205 


Uptakes (Smoke pipes) 


130.51 


206 


Propulsion Control Equip. 


10.97 


210 


Fuel Oil Service System 


10.10 


211 


Lubricating Oil System 


31.24 


250 


Propulsion Repair Parts 


8.50 


251 


Propulsion Operating Fluids 


42.23 




Total 


789.18 



Dry Weights of Principal Propulsion Components 



COMPONENTS NO. PER 


SHIP 


WT(L0NG TONS) 


Propulsion Gas Turbines 




4 


81.25 


Propulsion Reduction Gears, Inc. Acc 


* 

■ 


2 


149.56 


Propulsion Bed Plates 




2 


53-58 


Shafting *N/A(Total) 


156.43(34 ton/ft) 


CRP Propellers 


* 


2 


46.88 


CRP Propeller Hyd.Oil Power Module 


* 


2 


5.45 


CRP Propeller Oil Dist. Box 


* 


2 


2.32 


Line Shaft Bearings 


* 


5 


9.72 


Prop.GT Enclosure Cooling Fans 




4 


2.32 


IR Suppression Booster Pumps 




2 


1.61 


FO Service Booster Pumps 




4 


1.25 


L0 Service Pumps 




4 


4.46 


Prop GT Lube Oil Storage & 
Conditioning Assembly 




4 


2.95 


Prop GT Free Standing Electronics 
Enclosures 




2 


.38 


Total 






516.16 



♦Removed or replaced for electric propulsion 



Table 1.2 
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propulsion plant can be based between two ships identical 
in mission performance with differences only in propulsion 
plant size, weight and operating charcteristics. 

With this data as input, a ship synthesis model 

MM 

computer program was used to simulate each of the 
different propulsion systems and resultant ship. A 
sample output of this computer program is shown in 
Appendix H. , 

It soon became evident that the DD963 would not 
fit directly into the ship synthesis model without some 
"bias** being fed into the computer program. The synthesis 
model is based on all past design practices and limita- 
tions; the DD963 is based upon a present and new design 
philosophy. ^ 

The DD963 is a very spacious and roomy ship by any 
design standards. When the DD963 specifications were 
fed into the computer model, it generated a ship 900 long 
tons lighter than and 100,000 cubic feet of internal 
volume smaller than the actual DD963. This is due primarily 
to the DD963 being a weight limited ship with excess 
volume available for all shipboard functions. The "bias" 
had to be fed into the computer program to include this 
excess volume in the simulated ships. 

During the comparison of the electric propulsion 
system to the mechanical propulsion system, care had to 



28 



be taken not to lose this excess volume when the electric 
propulsion was inserted into the program. Loss of this 
excess volume would give a false indication of the 
desirablility of the electric propulsion systme due to 
a much smaller volume requirement for propulsion machinery. 

As a check to ensure that none of this excess volume 
was lost or misplaced by the computer program, the decision 
was made to make two comparisons of the electric propulsion 
system and the mechanical propulsion system. This will 
indicate if a volume saving is due to the electric pro- 
pulsion system or if it is due to the tightening up of 
a loose design. 

The first comparison is to be made between DD963S . 
with the excess volume as simulated in the computer with 
bias included in the program. The second comparison 
is to be made between DD963S as simulated on the computer 
based on past design limits where no excess volume is to 
be found. If the same relative volume and weight changes 
are observed in both comparisons, it would be a good 
indication that resultant weight and volume savings, (as 
simulated in the computer) , of the electric propulsion 
system based on the actual DD963 would indeed be a realized 
saving and not a false indication of volume reduction. 

The DD963 is limited in weight, with any weight decrease 
bringing about an improvement in overall ship characteristics. 
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The big advantage in the electric propulsion system is 
the weight savings due to decreased fuel and propulsion 
plant weight. The actual results of these comparisons 
can be found in Chapter 5» with sample computer outputs 
given in Appendix H and £. 
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CHAPTER 2 



PROPULSION MACHINERY DESCRIPTION 



2,1 Machinery Characteristics 



Once the prime mover and the transmission system 
have been selected the next design step is the selection 
of the basic propulsion plant configuration. At this 
point, the physical location of the individual pieces of 
equipment within the ship is not critical. The importance 
lies in how the power is to be passed from one unit to 



another to end up at the propeller with the least power 
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loss. A straight conversion of the DD 963 power plant to 
electric propulsion is shown in Fig. 2.1. 

The four gas turbines are "of type" LM-2500 configured 
in an enclosed module for marine installation. The 
enclosed module provides engine cooling, sound attenuation, 
internal lighting, view windows, and fire extinguishing 
capability. Output power is via two flexible couplings 
connected directly to the input shaft of the superconduct- 
ing generator. The characteristics of the gas turbine 
are given in Table 2.1. The performance characteristics 
are shown in Fig. 2.2 and 2.3* 



Rating Conditions 

Inlet Air Temperature (Power rating) 100°F 

Inlet Air Temperature (SFC rating ) 80°F 

Maximum Power Rating 



Brake Horsepower 
Power Turbine Speed 
Specific Fuel Consumption(SFC)Max 
Off Design Performance 



21.500 HP 
3,600 RPM 
0.42 lb/hp-hr 
see Fig. 2. 2 



Dimension (Module) 



Length 

Width 

Height 



26 ft. 6 in. 
9 ft. 0 in. 
9 ft. 6 in. 



DD 963 PROPULSION 
TURBINE CHARACTERISTICS 

Table 2.1 
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The power out of each gas turbine is transmitted 
directly to a 20,000 HP superconducting generator. The 
operating characteristics and dimensions of the generator 
are shown in Chapter 4. The output frequency of the 
generator will be 60 HZ when the turbine is at 3^00 RPM 
(full power) , and 30 HZ when the turbine is at 1800 RPM 
(min. operating speed of power turbine), see Fig. 2.3* At 
turbine speeds between I800-3600 RPM, this frequency 
change will directly control the speed of the motor and 
hence the propeller speed. 

The switching power control unit produces one of the 
big advantages of an electric propulsion system over a 
conventional mechanical system. In the mechanical system, 
any given gas turbine is physically connected to a specific 
propeller. In the electric system, any gas turbine-genera- 
tor can be connected to either propeller. At low speeds, 
this can result in great fuel savings, as both propellers 
can be driven by one gas turbine. A single gas turbine 
operating at higher power is much more fuel economical 
than two turbines operating at a lower power rating. In 
Fig. 2.2, one turbine providing 17,000 HP has an SFC of 
.445 lb/hp-hr at 3000 RPM, while two turbines providing 
17.000 HP (8, 500 HP each) have an SFC of .56 lb/hp-hr at 
2300 RPM. This works out to be a savings of 1955 lb/hr. 
This is only one example of savings offered by an electric 
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SHAFT HORSEPOWER 



Figure 2 . 2 - ESTIMATED MINIMUM PERFORMANCE 
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SPECIFIC FUEL CONSUMPTION lb/hp-hr 







I . — ..I. —»———— I I .. I t I 

o 3000 iqooo 15,000 2 0,000 25,000 

SHAFT HORSEPOWER 



SPECIFIC FUEL CONSUMPTION OF THE 
LM-2500 GAS TURBINE ENGINE AS A 
FUNCTION OF POWER LEVEL AND 
OPTIMUM POWER TURBINE SPEED 



Figure 2.3 
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propulsion system. A detailed breakdown of overall 
savings are given in Chapter 6. 

The control unit performs several important functions. 
The most important is providing the necessary reversibility 
of the propeller by controlling the direction of rotation 
of the superconducting motor. In addition to controlling 
the direction of the motor, the control unit can also 
control the speed of the motor. Speed control is accompl- 
ished by controlling the frequency of the electrical power 
to the motor. A possible method of propeller control 
for propeller speeds of 85 to 170 RPMs, would be the 
control unit employing a fixed frequency reduction and 
the speed of the propeller being controlled directly by 

4 J 

the speed of the gas turbine (I8OO-36OO RPM) and the 
electrical output of the generator (30 to 60 HZ) . At 
propeller speeds of 30 to 85 RPMs, the control unit 
would control the speed of the motor by using a variable 
frequency reduction. For this range of propeller RPM, 
the gas turbine would be held at a constant speed with 
the generator delivering power at a constant frequency. 

In this example, the control unit controls both the speed 
and direction of the propeller. Examples of different 
control units have been proposed. ( 2 )(3)(5)(7) 

In basic, design, there is no difference between 
a motor and a generator. Superconducting motors are 
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very similar to the superconducting generators. The major 
difference being, the motors are 40,000 HP each at 200 RPM 
maximum design speed, where the generators are only 

20,000 HP each at 3^00 RPM maximum design speed. 

Last in propulsion is the propeller which delivers 
the output of the motor to the water. The controllable- 
reversible-pitch propeller presently on the DD963 can 
be replaced by fixed pitch propellers, which provides a 
higher effeciency. See Table 2.2. 



Rated Power 
Rating RPM 

Effeciency (open water) 

Diameter 

No. of Blades 

Hub Ratio 

Expanded Area Ratio 
Weight 



CRP 

40,000 HP 
168 RPM 

7 Ofo 

17 ft. 

5 

.30 

.75 

23.4 L Tons 



FIXED PITCH 

40,000 HP 
170 RPM 

73 % 

17 ft. 

5 

.17 

.75 

18 L Tons 



PROPELLER CHARACTERISTICS 



Table 2.2 



The cryogenic refridgeration system is not shown in 
Fig. 2.1 , but it is an essential part of the propulsion 
system. The best configuration at present is to have one 
refridgeration system per electric motor, and at least 
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one system per two electric generators. With crosst- 
connect piping installed, this will provide the required 
reliability without excessive redundancy in the cryogenic 
systems . 

Chapter 3 contains the design of the superconducting 
motors and generators described in the above propulsion 
system. The remainder of the electrical transmission 
system is covered in Chapter 4. 
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CHAPTER 3 

DESIGN OF SUPERCONDUCTING MACHINE 
3.1 Introduction 

An optimization design program for the basic electri- 
cal design of superconducting generators and motors, has 
been used in this thesis, see Appendix A for computer 
program listing. This computer program optimizes the 
machine by seeking a design which minimizes a total "cost 
function", which is a function of machine weight, dimensions 
and operating characteristics. The operation of the 
optimization program and the bulk of the subroutine CF 
can be found in KIRTLEY et. al.^ 2 ^^ The portion of 
the subroutine CF which computes the damper stress and 

/ON 

thickness is from a thesis by Furuyama. ' A detailed 
development of the mathematical equations and theory 
abstracted here, can be found in these previous works. 

The function of the computer program is to select 
an optimum design of a superconducting generator or motor 
for a machine of a given physical configuration and horse- 
power rating. An initial set of dimensions is assumed 
for the machine. The computer program, using this 
initial guess as a starting point, attempts to design 
a machine that provides the required horsepower, with 
the least weight, smallest dimensions and least internal 
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power loss based on a set of rules. 

The program is based on an optimization approach. 

For each trial design, a single number, called a "penali- 
zed cost" is generated. The program attempts to find a 
design for which this "penalized cost" function is a 
minimum. To accomplish this, it searches over the values 
of seven machine dimensions and the value of field current 
see Table J.l. The cost function is the product of two 
numbers : 

1. The "cost" is the sum of material weights. 

2. Penalty functions are established for 
several variables that will have values 
that are either acceptable or not acceptable 
The penalty functions have very large values 
when their associated variables are unaccept 
able . The penalty functions are multiplied 
together and then multiplied by the "cost" 
to form the "penalized cost" . Penalty 
functions exist for: 

Maximum field current density 

Shaft critical speed 

Shaft stress 

Shield Flux limit 

Damper stress 

Armature inner radius 

Armature insulation thickness 
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EXIT 



BLOCK DIAGRAM OF PROGRAM 
Figure 3.1 
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See eq. 3*91 for the generation of the penalty functions. 

The optimization approach works by taking each of 
the eight search variables separately, and attempting to 
find a local minimum of the value returned by CF for 
each variable. Three calls to CF are made for each 
variable, with the value of the variable incremented 
twice by a fixed value: 



Y 1 = CF(V q ) 

Y 2 = CF(V Q + T^ DV) 

Y_ = CF(V + 2 x T x DV) 
j o w 



(3.1) 

(3.2) 
(3-3) 



A second order curve of the following form, when fit to 
these three points 



Y(V) = aV^ + bV + C 



(3-*0 



will have as its postition of zero slope 



* y*i - 4y 2 + y 3 

V = V o + T w DV 2Y, - 4 y! T 2 t 



(3-5) 



This resulting position will be a minimum of that second 
order curve if the second derivative of Y with respect 
to V is greater than zero 



The optimization routine then selects the "optimum" 
machine design based on the penalty functions generated 

I in subroutine CF, which is fit to the curve of eq. (2.4) 

to find the minimum "cost" times "penalty". Following 

is a summary of the calculations performed in the main 

subroutine CF. 

3.2 Subroutine CF 
— 

Subroutine CF utilizes two subroutines, CS and CM, 

which calculate the geometric parameters used in the 

S 

inductance expressions. 

The optimization variables or search variables, (see 
Table 3.1) 



R fi 


is 


field inner radius 


T hf 


is 


field thickness 


G fk 


is 


field- to- damper gap 


T hk 


is 


damper thickness 


G ka 


is 


damper-to-armature gap 


T ha 


is 


armature thickness 


G as 


is 


armature-to-core gap 


X f 


is 


field current density 



Table 3.1 

are passed into CF through the vector V, which is the 
only argument into CF from the optimization program. 
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With the exception of the values of V and the returned 
value of CF, which is the penalized cost CF, all other 
variables used by CF are fully self-contained within 
the subroutine. All other values generated by CF and 
passed to the output are the optimum machine dimensions 
and characteristics. 

3.2.1 Machine Length and Synchronous Reactance 

The calculation of machine length and synchronous 
reactance is generalized by an arbitrary number of armature 
phases and the inclusion of an armature winding factor. 

This calculation is complicated because of a voltage drop 
in synchronous reactance over the unknown machine length, 
and by the effect of the end turns. The machine rating 
in volt-amperes is given by: 



VA 



N, 



i*a V t ] 



( 3 . 7 ) 



where 



VA 

•Va 



is rating in volt-amperes 

is the number of armature phases 

is phase voltage (RMS) 

is phase current (RMS) 



The RMS value of the internal voltage may be written as: 



E 



f 







( 3 . 8 ) 



where 




<i> 



is internal voltage 

is field-armature mutual inductance, 

given by eq. ( 3 * 13 ) 

is field current 

is electrical angular frequency 



The relationship of to can be obtained from the 
phasor diagram for machine operation (see Fig. 3*1)* 




for 

Operation as a Generator 

) 

Figure 3.2 

The law of cosines is applied to the above diagram to 
yield: 

E f 2 = V t 2 + x A 2 I t 2 + ZV t x A I t sin* (3-9) 
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where is machine synchronous reactance in ohms; 





(3.10) 



and t is power factor angle. L is the phase self -induct- 

H. 

ance, given by eq. (3.14). By dividing and then rearrang- 

p 

ing eq.(3.9) by E f , the ratio between terminal voltage 
and internal voltage is obtained: 



(J 6 ) = g 

E f ^ 



2 2 

1 - x„ cos f - x„ sin ^ 
a a T 



where 



x_ = 



X a I t 



a E, 



(3-11) 



(3.12) 



is the synchronous inductance normalized to internal 
voltage. The two inductances used here are: 



H = 



_ 32 * a » 0 H a M f sin(^p) sin(-^p)(l- i i p+2 )[R 



P "Ve e wfe (1 -y > (1 - x > 



'IT 




C k 
m wa 



(3-13) 



I 



L = 
a 



164 



u N 

oa^o a 



„2, wae\ 
sm l — 5 — ; 



P «Ve 2 f 1 -* 2 ’ 2 



c k * 

s wa 



(3*1*0 
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where 



*a 


is 


^oa 


is 


. 0 


is 


N a 


is 


N f 


is 


0 

wae 


is 


9 wfe 


is 


X 


is 


y 


is 


R fo 


is 


R ao 


is 


C m 


is 


P 


is 


C s 


is 


R wa 


is 



active length for mutual coupling 
active length for self- inductance 
permeability of free space 
number of armature turns 
number of field winding turns 
armature phase winding angle 
field winding angle 
armature radius ratio 

field radius ratio 

field outer radius 
armature outer radius 
the mutual coupling coefficient see eq.( 3 * 30 ) 
the number of Pole pairs 

the self- inductance coefficient see eq.( 3 « 31 ) 
the armature winding factor 



x = 



y = 



ai 



ao 

l fi 



'•fo 



The field and armature currents are related to current 
densities by: 



h - 



J f 0 wfe R f 0 2 (1 -y 2) 



2N 



ft 



I a = 



J a Ve 8 ao 2 ^ 



2N. 



at 



( 3 . 15 ) 



( 3 . 16 ) 
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i 



i 



l 



Now, substituting eqs. (3-8) , (3-H) . (3.13) . (-3.14) , (3*15) 



and (3.16) into eq.(3.7) and performing some manipula- 
tions yieldsi 



VA = P. 



p a 



2 2 

l-x„ cos 0 - sin^ 
a T a 



(3.17) 



where 



p p = 



8 % a“ -0 V f q-y p, ' 2 )R f0 2+p R - 2 " p C 



ao 



\1 2-tt 



k 



wa 



ft 



sin(^fS) 8in(^fS) 



(3- 18) 



Substituting eqs. (3*8) , (3*10) , (3*13) 1 and (3.14) into 
eq. (3*12) yields: 



x a =x il 



oa 



(3-19) 



where 



| 



x i 



£ %a J a s^^r> 

4 



32 + p 



J f sln( 9 wfe' 



ao 



L R foj 



s wa 



C m (l-y 2+P ) 



( 3 . 20 ) 



The effect of end turns on self- inductance can be 
represented by: 



t - 1, + 4 

oa a 



(3.21) 



, 
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The end turn correction length a is postulated to be: 



A = 



R + R . 
ao ai - 



R fo + R fi 



c b£j 



L bX 



( 3 - 22 ) 



In this equation is used to describe different end 
turn forms and the factor use< * assign part of 

the end region to mutual coupling. 

Combining eqs. (3*19) , (3*20) and (3.21) yields: 









a +a 


2 


& + A 


sin^ 


VA = P v ( & 




!-Xi 2 


a a 


2 , 

cos i - x^ 


a 


\ 


L*aJ 


L«aJ 









(3.23) 


By defining: 


« - VA 

2 P P 


(3.24) 


and then dividing eq.(3«23) by VA, the following 


is 


obtained: 






1 =Ja 2 


- x^ 2 (a+P) 2 cos 2 ^ - x^(a+£)sin^ 


(3.25) 


where 


a - 7T" 
*z 


(3-26) 




a _ A 

e “ T 

A z 


(3-27) 
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Eq.(3*25) can be solved for a: 



a 



a + N a 2 + 1 + 2x^Psin^ + x^ 2 p 2 



0 . 28 ) 



where 



a 




(3.29) 



Per-unit synclironous reactance based on terminal voltage 
is then determined to be: 



3.2.2 Geometric Coefficient C and C 
m s 

The geometric coefficient is calculated by the 
subroutine CM, with arguments p,x and w. For mutual 
inductances involving the armature winding, C m is found 
as follows: 



x 



d 



E 



f 




C m = -g(-ln x + 4 ( l — x 2 ) w^ ) if p = 2 



C 



m 



1 ~x 2 " P + - j^ (l-x 2+ P)w 2 P 

4-p 2 



if p/2 (3-30) 
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The geometric coefficient C g is calculated by CS with 
the arguments p,x and w for self inductance. 






if p = 2 



° s = 



(2-p)4x pf2 + 3x p + 2(|^)(1 -x 2+P ) 2 2p 
p(4-p 2 ) 



if p / 2 

(3.3D 



The above procedure is used for the self- inductance 

| 

of all windings, by substituting appropriate parameters 
for x and w. 

Overall lengths for the damper and field winding 
Eire computed in the same fashion as for the armature. 
Bearing length is assumed to be damper length plus the 
length of the thermal distance pieces. The thermal 
distance piece length LTH, is an input constant. 

1 3.2.3 Effective Current Density 

! 

■ As stated, the armature conductors are not aligned 

with the axis of the machine; therefore, the axial com- 
ponent of the current which produces the interaction is 
found by: 



J ah = J a cos e h 



(3.32) 



51 



where J is the total current density in the helical 

cl 

path, and is the helix angle: 



©h = tan' 



-1 " R ai 



( 3 . 33 ) 



If x^ is computed according to eq.(3.20) and using 
total current density J , the internally based synchronous 

cl 

reactance will be found to be: 



x = x. cos©, (1 + — ) 
a l h 

a 



where a is given by eq.(3*22). 

The machine length x^ is found by: 



L = 



YA 



Pp COS0 h ( E ^) 



( 3 - 34 ) 



( 3 . 35 ) 



where is found by eq.(3«18) and using total current 
density J . 

a 

3.2.4 Transient and Subtransient Electrical Parameters 
Transient reactance for the machine can be found 
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by: 



( 3 . 36 ) 



I 



d 



1 - 



M 



L a L f. 



The equation for L^. is the same as eq.(3.l4) for L a with 
the substitution of y for x, R fo for R &o and 9 wfe for 
9 . Substituting these values into eq.(3.36) yields: 

W3.6 



x d = x d 



1_4 a wa 



C m 2 ( l " yPf2)2 



^oa^of C s( p,x,R ao/R ) C s( p,y,R fo/fc ) 



R fo 2p 
ao 



( 3 - 37 ) 



Subtransient reactance is given by: 



x d" = x d 



1-2 ( 



i 2 C 2 

a .) 111 



4a*of c s <P.x.w) C s( p, Z ,R koA j 



(^ S ) 2 P d -^ 2 ) 2 

R ao 



( 3 . 38 ) 



The dynamic performance of superconducting machines is 
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I 

i 







a very important aspect of the overall machine design. 
This is brought about by a conflict between adequate 
rotor shielding, which requires a high rotor conductivity 
and a damping of rotor swings, which requires a lower 
conductivity. The damper time constant T g and the 
armature time constant T are given by; 

T s - f '‘o R ko (R ko - R ki> °k (1 + R kiA k0 )2P (3 - 39) 



T = 



2'o sin 2 <%^) V a R ao 2 C S <P’*- R ao/R “wa" 

s 



(3.40) 



The open-circuit subtransient time constant is: 



rp «• 

1 do 





(3.41) 



3.2,5 Field Current Rise 

An estimate must be made of transient field current 
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resulting after a fault or short circuit. The value 
calculated here, corresponding to maximum field current 
during a critical swing, is used later to determine if 
the field current is within limits. 

Two assumptions are made to simplify the calculations: 

1. Field flux is constant over the period of 
the swing. 

2. There is no shielding of the field for 
this transient. 

is voltage behind synchronous reactance corres- 
ponding to the operating condition in which: 

I^ o is field current 

6 q is torque angle 

1^1 and 6^ correspond to the fault condition 
The field current rise is: 



I fl _ _2_ x d ~ x d . , , 

X fo ~ E fo V + x e ° 



(3.^2) 



, The initial torque angle is found by: 



6 0 = sin' 



-1 



p(x d + x e ) 



V E 



fo 



( 3 .^ 3 ) 
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2 

The values of E^ o and are f ound from the following 
expressions: 



E fo 2 = V t 2 + ( x d J t )2 + 2 Wt sin * 



( 3 - 44 ) 



V. 



2 



= V 



2 




( 3 .^ 5 ) 



00 



t 



3.2.6 Torque Tube Thickness 

The inner member of the rotor is a heavy- walled 
cylinder called the torque tube. The primary requirement 
on the torque tube is that the stresses be kept below 
the yield stress of the material. This yield criterion 
is necessary in order to avoid any fatigue failures or 
changes in dimensions which could result in a mechanic- 
ally unbalanced rotor. 

The computer program attempts to estimate the 
thinnest torque tube that will take the worst-case 
fault torque duty. During a fault, the radial forces 
generated will force the torque tube into an out of 
round condition. All circumferentially dependent loads 
are averaged and treated as uniformly distributed forces. 

The max shear stress is given by: 



T. 



max 




( 3 . 46 ) 



56 



The two components of stress are shear stress resulting 
from torque 



r =■ 
S 



tt(R 



2Tr 

"T 

ho 



- O 



0.47) 



and tensile centrifugal stress 



3+ v 2 
°6 = 8 “ 9 “ 



Rx.. 2 + Rx.. 2 + R hp 2 R hj 2 - l+3 v _2 



bo 



bi 



3+v 



(3.48) 



where r is the radius 

is 

U± 

R 



R bi 



bo 






inner radius of torque tube 
is outer radius of torque tube 
p is mass density for torque tube 

v is Poisson's ratio for torque tube 

T is torque 

The computer program tests the maximum stress 
eq.(3.46) at the inner and outer radii of the torque 
tube. It compares the larger of the two stresses with 
the stress limits which are an input to the program. The 
program will accept torque tube stresses between: 
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.95 



0 . 49 ) 



T T r 

lim max lim 



I 

I 

i 



If the maximum stress is not within these limits, the 
program computes a new value of inner torque tube radius 
divided by the outer torque tube radius which is a ratio 
of torque tube thickness and then recomputes the maximum 
stress limits. 

If eq.(3.^7) is not satisfied after 250 tries, the 

programm assumes a solid shaft and computes the maximum 

stress in this shaft. The maximum stress level r 

max 

will later be used in a penalty function for torque tube 
stresses . 

3.2.7 Damper 

The damper is a thin conducting cylinder located 
at the outermost diameter of the rotor. The principal 
functions of this damper shield axe to shield the super- 
conducting field winding from alternating magnetic 
fields and to damp the mechanical oscillation of the 
rotor. In the event of a terminal fault, the damper 
has to withstand the strong crushing and torque loads. 
These loads must be computed in order to insure adequate 
thickness of the damper. The force per unit area on the 
damper during a terminal fault is given by: 



(3-50) 



a 



r 





) 



a e = 



'o< H eo 



- H ei> H r 



(3.5D 



In a magnetic field, where H r is the radial component 
and the tangential components inside and outside the 

/ON 

cylinder are H 0 . and H- respectively. 1 

Ul U O 

The magnitude of the fundamental component of the 
tangential field at the damper due to the average rms 
armature current density J is given by: 

cl 



KL = 



lOT 

TT 



sin(- 



wae 



> R ao 1-^ + 3 <X-x j )( r 



R a „ 2l 

mi J 



(3.52) 



The tangential field at the angle 6 outside the damper 
due to the armature current immediately aftep a fault is: 



H ea = ^ Sin <-F> R ac J ar 



-1 o 2 

1-x + i(l-x 3 )(=a2) 

3 mi J 



V, _ 

u |j;os(wt+p(-e)-cos(^-e)j+ cos(wt+pf-e+^- ifO 



L A d 






(3-53) 
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The tangential field inside the damper due to the rated 
armature current, which is not affected by the fault, 
is given bys 



H 



eao 



= iJJ 



TT 



sin(- 



e 



wae 



>R ao J ar 



1-x + ±(l-x- , )(# 2 ) 



mi. 



cos(wt + 0'-e+-^-f) 



(3.54) 



The tangential field at the damper, due to the field 
current, which is constant before and immediatley after 
the fault is: 



"of " ji sinf^R^Ml-y 3 ) 

s 

cos(wt + 0* - 0 + Cf) 



R c 2 
1 * ( 5 a 7 ) 

mi _ 



’f 

(3-55) 



The current in the damper is induced so that the radial 
component of the total field is kept constant before and 
after the fault. Immediately after the fault, the field 
due to this induced current is: 



H es = 



R„ 2 
l-(tr-) 
R mi 



R 2 
l+Ur 8 -) 

l_ R mi _i 



# sin(^)R ao 



1 -5 R -„ 2 

1-x + i(l-i 3 )(jS2) 

3 R mi J 



J V t 



ar x d ,n 



[cos(wt + 0 - e) - costp' - efj (3.56) 
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Adding up these components, the total tangential 
field outside the damper immediatley after a fault is 
obtained. 



H a = + H a + H 

eo ea es ef 



The total tangential field for the inside is: 



H ei = H ef + H eao 



( 3 . 57 ) 



( 3 . 58 ) 



To obtain the maximum value of a r , it is assumed that the 
maximum a r occurs when the traveling wave (sum of three 

terms which have «>t in eq.(3*57)) comes to the same phase 

( 8 ) 

as the standing wave (cos(0-e) terms in eq.(3»57). 1 



o T = cr rl + a r2 cos 2(0 - e + y) (3*59) 



where 



d rl = ’T (f2 ' d2 ’ e2) 



*r2 T \ 



(F 2 - E 2 + D 2 ) + 4E 2 D 2 



_ 1 , -1 2ED 

Y p tan ? 2 2 

* F - E + 



(3.60) 



( 3 . 61 ) 
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A = H ^ sin 4 , + 
ao Y 



2H 



a 



1 + ( R § “ ) 
mi 



~2 + H f cos 6 



(3.62) 



B = H ao cos^ + H f sin 6 



(3.63) 



2H 



C = 



a 



R a 2 

1 + ( R^" ) 
mi 



(3.64) 



D = H ao cos(tan _ ' L j + <p ) - H^sin(tan -1 ^ - 6) 



(3*65) 



E = -H„^ sin(tan _ '*‘ ^ + <p ) - H^cos( tan - "* - ^ - 6) 



ao 



= 4 



A 2 + B 2 - C 



( 3 . 66 ) 

(3.67) 



M = —■ (COS a - -) 
c. TT 



( 3 . 68 ) 



For concentrated radial force as shown in Fig. 3«3» 
the radial displacement u and the bending moment M are 
given by: 




M_ 

El 



(3.69) 



where M is found from eq.(3.68). For the bending stresses: 



6 _ M _ PR 

b ~ Z " 2Z 



(cos or 




(3.70) 
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Concentrated Radial Force 
On The Damper 

Figure 3.3 




cos 2(p-a) 



Distributed Radial Force 
On The Damper 



Figure 3.4 
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The deflection is then solved: 



Pr 3 pr 3 _ p R 3 
u “ ttEI “ Vei a sinrt “ Wl cos “ 



(3.7D 



For the distributed load as shown in Fig. 3 . 4 , the 
deflection at angle « is: 

,4 



4 ^ 

u(«) = -s ~ — cos 2 « 

3 E t 3 



(3.72) 



where t is the thickness of the damper. 

The bending stress in Fig. 3*4 is then given by: 



2a r2 R 

a, = s- COS 2 a 

b t 2 



( 3 . 73 ) 



Adding the deflection and the stress due to the centri- 
fugal force and the uniform magnetic force to eqs.( 3 * 72 ) 
and(3.73) the total deflection and stress are obtained. 
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2 cr 



r2 



total 



cos 2« + p " 2 



(3-74) 
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(3.75) 
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Using eq. (3*7*0. "the thickness t of the damper can 
be calculated such that the maximum allowable stress and 
deflection will not be exceeded during a terminal fault. 



3.2.8 Negative Sequence Losses 

Negative sequence currents in the stator produce a 
magnetic field distribution which rotates in a sense 
opposite that of the rotor. Since the damper is a good 
shield, it excludes this magnetic field from the rotor 
and must match this field at its surface. The magnetic 
field at the surface of the rotor is given by: 



The current density in the damper is: 




The negative sequence loss in the damper is then given 




(3.76) 



P 



sh 




(3-77) 
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where 






(3-78) 



3.2.9 Armature Losses 

Conduction losses in the armature are calculated by; 



P a - 0 ± 0 wae % R ao 2 C 

3 . 3 



( 3 . 80 ) 



3.2.10 Field at Shield Radius 

This is the radial magnetic flux density at the 
inner radius, R , of the stator core: 

S 



wfei 



B 



J f sin(-=±*)_ R fo .2+p 



rs 



TT 



( 2+p) 



— R s ( ^) ( 3-8D 



The core outer radius is then found by using the 
magnetic flux density and the inner radius: 



R 0 „ = R (1 + j: 
SO S B 



rs 



smax x 



(3.82) 



3.2.11 Field at an Inner Corner of Field Winding 

The highest field intensity is assumed to occur at 
the inner radius of the field winding, and at positions 
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adjacent to pole faces. The radial and azimuthal fields 
are calculated by: 



H r ' § H n W <3- 8 3) 

<3.84) 



where 



2J„ 2+np « 

H n = IT — R fi y p c m (n P’y' R foA ) <3-85) 



The maximum field intensity is then: 



*W =4«r 2 + H 0 2 



( 3 . 86 ) 



3.2.12 Rotor Critical Speed 

The rotor is assumed to be a beam of constant stiff- 
ness and mass, simply supported at both ends for the 
purpose of estimating rotor first-critical speed. The 
first critical frequency is then: 



• c = 9.875 



N 



El 

1 W 

*br 



M 



br 

M 



IS 

is 



bearing length 
mass per unit length 



(3.87) 
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It is further assumed that the only stiffness is 



provided by the torque tube: 





(3-88) 



The mass per unit length includes the entire rotor 
which consist of the torque tube, field windings, hoop 
binding material for field winding, shield and damper. 
The length is assumed to be the bearing length • 
3.2.13 Stator Core Losses 

It has been assumed that all of the core is operat- 
ing at roughly a uniform flux density, and hence a uni- 
form loss density. 



P = M P 
core core m 



where : 

M 

core is the total core mass 

P^ is dissipation per unit mass 



(3-89) 



Core dissipation per unit mass is estimated to be a 
simple power function of flux density: 



B V 

p„ = f 0 i r ss -'> 



m 



smax 



where : 

B is computed in eq.(3«81) 

rs 



(3-90) 
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B „„„ is core flux limit 

smax 

P Q is dissipation per unit mass when the 

core is operating at its limiting 
flux density 

Y reflects the rate of change of loss 

with flux density 

All the above, except for B , are input variables. 

r S 

3.2.14 Cost Function 

The non-penalized cost is the cost of the machine 
in weight. The weight includes the support tube, damper, 
armature windings, field winding, iron in the core and 
binding material. A weight loss in KG/watt is also found. 

3.2.15 Penalty Functions 

The cost of the machine, such as weight, current 

i 

losses or stresses is modified by multiplying by a set 
of penalty functions. These are of the form: 

F p = .9 + .KQ/Qjj ) 15 (3-91) 

where : 

Q is quantity being penalized 

Qj is the maximum limit for that quantity 

This quantity is close to 1 for values of Q less than 
Q| , but becomes very large for values of Q greater than 
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The field- current-limit penalty function is slightly 
different. The extremes of magnetic field and current 
density in the field winding are H maY eq.(3»86) and J.j.1^ 
eq.(3*^2). This combination must be compared to the 
H-J curve. The magnetic field-current density curve 
is approximated by a six-segment piece-wise-linear curve 
as shown in Fig. 3 * 5 * The data representing this curve 
is input at two five-element vectors, one representing 
the values of H, and the other representing the vlues 
of J. 



H-J CURVE 




If H is greater than H, , the field- current-limit 
max jl 

penalty function is set to a very large number. If H 

max 
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falls between the values of and H^, the computer 
program does a direct linear interpolation to find the 
critical current density J. If H wax is less than , 
the critical current density is set to . 

The final penalized cost function CF, is the product 
of the cost and all of the penalty functions. This is 
the value returned by CF to the optimization main program. 
The value of CF is used as a figure of merit for each 
different machine iteration. The machine selected by the 
program is the one with the lowest value of CF . 



CHAPTER 4 



Superconducting Propulsion Plant 
4.1 Introduction 

Chater 1 has covered the basic machinery configuration 
of the electric propulsion system. The actual design 
of the superconducting machine is described in Chapter 3. 

The "optimized" electric machines and the associated sub- 
systems necessary to complete the entire propulsion 
system will be discussed here. 

The output of the optimization program is an "optimum 
machine" only for one particular set of assumptions as to 
machine design requirements and costs. The fixed inputs 
common to both the generator and the motor are given in 
Appendix F. The initial guess for the search variables 
or optimization variables are shown in Appendix D. The 
combination of these inputs yield the optimized generator 
design and motor design described in Sections 4.2 and 4.3 
of this chapter. The requirements of the machines described 
in sections 4.2 and 4.3 are then used to determine the 
characteri sties of the remaining subsystems in Section 
j 4.4. 

4,2 Superconducting Generator 
' The first step of the electric propulsion system 

design process was the design of a 20,000 HP 3,600 RPM 
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synchronous generator driven by 20,000 HP gas turbines. 

A summary of the results sure shown in Table 4.1. The 
detailed computer output results are contained in 
Appendix B. 

All of the machine dimensions, machine volume and 
the weight of active parts are computed by the computer 
program. A steel outer shell for machine structure 
support was included as additional weight and volume, 
resulting in the total machine volume and weight. (The 
computer program did not compute the outer shell character- 
istics. Provisions can be made at a later date to include 
this in the program.) The steel outer shell was assumed 
to be one-half inch thick (.012? M) plus the structural 
supports required to mount the machine to a foundation. 

The end bells were assumed to be j/b inch thick steel 
and the same diameter as the outer diameter of the newly 
added steel shell. ' 

Using these new dimensions, the new volume is calculated 
to be 1.1 B/P and the new total weight is 4309.2 KG. The 
weight of 4309.2 KG also takes into account the weight 
of the generator feet and steel supports necessary to 
mount the machine to its foundation. 

4.3 Superconducting Motor 

The next output of the optimization design procedure 
was a 40,000 horsepower, 200 RPM synchronous motor. A 
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A SUMMARY OF RESULTS 
FOR 

A 20,000 HP GENERATOR 
Table 4.1 



Mechanical Rating 


20,107 HP 




Electrical Rating 


15 MVA 




Mechanical Speed 


3.600 RPM 




Number of Poles 


2 




Active Length 


.69 M 


27.2 IN 


Overall Length 


1.52 M 


59.9 IN 


Field Winding Inside Diameter 


.19 M 


7.7 IN 


Outside Diameter 


.25 M 


9.8 IN 


Armature Winding Inside Diameter 


.43 M 


16.9 IN 


Outside Diameter 


.62 M 


24.4 IN 


Iron Shield Inside Diameter 


.67 M 


26.4 IN 


Outside Diameter 


.92 M 


36.2 IN 


Machine Volume 


1.01 M 3 


35.7 FT 


Weight of Active Parts 


3,175.4 KG 


3.1 T 


Total Machine Volume 


1.1 M 3 


38.8 FT 


Shell .EndBells, Bearings & 
Structural Support 


1,133.8 KG 


.9 T 


Total Weight of Machine 


4,309.2 KG 


4.0 T 


Synchronous Reactance 


1.38 




Transient Reactance 


1.18 




Subtransient Reactance 


1.20 
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summary of the results for the 40,000 HP motor are shown 
in Table 4.2. The detailed computer output for the motor 
is contained in Appendix G. 

The dimensions of the active machine parts sure 
computed by the optimization program, while the weight 
and volume of the steel structural shell and the weight 
of the motor feet were added to the computer output to 
produce the total machine weight and volume. For the 
motor, the steel outer shell was assumed to be 3/4 inches 
thick and the end bells 1 inch thick. The structural 
support of the motor is heavier than that for the generators 
because the motor must support the high torque loads 
associated with the propeller providing thrust to move 
the ship through the water. 

The summary of results for the 30,000 HP motor are 
shown in Table 4.3. The 30,000 HP motor was designed 
when it became necessary to reduce the size of the ship's 
propulsion plant from 80,000 shaft horsepower to 60,000 

i 

shaft horsepower. Chapter 5 Sections 5«3»1 and 5»3*2 
explain why the installed shaft horsepower was reduced. 

The total machine weight and volume were computed in the 
i same manner as that for the 40,000 HP motor. The computer 
output for the 30 » 000 HP motor is contained in Appendix D. 
4.4 Subsystems of Superconducting Machinery 

Cryogenic refrigeration systems are necessary to 
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A SUMMARY OF RESULTS 



FOR 

A 40,000 HP MOTOR 
Table 4.2 



Mechanical Rating 


40,214 HP 




Electrical Rating 


30 MVA 




Mechanical Speed 


200 RPM 




Number of Poles 


6 




Active length 


1.51 M 


59.4 IN 


Overall Length 


2.35 M 


92.4 IN 


Field Winding Inside Diameter 


.86 M 


34.0 IN 


Outside Diameter 


.93 M 


36.4 IN 


Armature Winding Inside Diameter 


1.11 M 


43.7 IN 


Outside Diameter 


1.46 M 


57.^ IN 


Iron Shield Inside Diameter 


1.52 M 


59.9 IN 


Outside Diameter 


1.73 M 


68.2 IN 


Machine Volume 


5.52 M 3 


195.3 FT ; 


Weight of Active Parts 


13.394.13 kg 


13.2 T 


Total Machine Volume 


6.01 M 3 


212.8 FT- 


Shell .EndBells , Bearings & 
Structural Support 


1948.17 KG 


1.9 T 


Total Weight of Machine 


15,3^2.3 KG 


15.1 T 


Synchronous Reactance 


1.01 




Transient Reactance 


.87 




Subtransient Reactance 


.86 





76 



A SUMMARY OF RESULTS 
FOR 

A 30,000 HP MOTOR 
Table 4,3 



Mechanical Rating 




30,161 HP 




Electrical Rating 




23 MVA 




Mechanical Speed 






200 






Number of Poles 






6 






Active Length 






1.19 


M 


46.85 IN 


Overall Length 






2.01 


M 


79.1 IN 


Field Winding 


Inside Diameter 


• 79 


M 


30.94 IN 




Outside 


Diameter 


.86 


M 


33.86 IN 


Armature Winding 


Inside Diameter 


1.05 


M 


41.15 IN 




Outside 


Diameter 


1.40 


M 


55.1 IN 


Iron Shield 


Inside Diameter 


1.46 


M 


57.56 IN 




Outside 


Diameter 


1.67 


M 


65.71 IN 


Machine Volume 






4.68 


M 3 


165.6 FT- 


Weight of Active 


Parts 


10 


,577-7 


KG 


10.4 T 


Total Machine Volume 




4.68 


M 3 


165.6 FT 3 



Shell ,EndBells, Bearings & 
Structural Support 

Total Weight of Machine 

Synchronous Reactance 

Transient Reactance 

Subtransient Reactance 



1,689.2 KG 
12,266.9 KG 
.98 
.87 
.86 



1.6 T 
12 T 
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provide the helium for super-cooling the previously 
designed motors and generators. A diagram of a sample 
refrigerator/liquefier system is shown in Figure 4.1. 

The liquefier produces liquid helium which is stored in 
the liquid accumulator. The liquid helium is piped 
directly from the liquid accumulator to the superconducting 
field winding where it provides the necessary cooling by 
expanding back into a vapor. The returning helium vapor 
from the field winding passes through a precooler in the 
Liquefier, increasing the overall cycle effeciency. At 
present, most of the high capacity helium refrigeration 
units in use have been constructed for fixed installations, 
where weight and size are of secondary importance. The 
fixed installation type of compressors are a low-speed 
reciprocating compressor of bulky construction. They 
generally require massive foundations because of the 
reciprocating loads. Machinery of this size is totally 
unacceptable for shipboard use, where weight and size are 
at a premium. 

The rotary compressor is a new entry to the field 
of helium refrigeration systems which will overcome the 
disadvantages of the reciprocating compressors. At 
present, there is not a great deal of information concern- 
ing high capacity rotary units. Experiments at MIT and 
elsewhere ^ ^ have produced similar results as to the size, 
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STATION 

1 

2 

3 

4 




HELIUM LIQUEFIER/frBFRIGERATOR DIAGRAM 
Figure *V.l 
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weight and operating characteristics of the rotary com- 
pressor necessary to supply the cooling requirements of 
the motors and generators designed in this thesis. 

The cooling requirement for the superconducting 
machines is an input to the program. This input, RP, is 
a penalty function for cryogenic refrigeration based on 
specific power consuption (watts input per watt at 4.2° K) . 
Based on experience at MIT, the value of RP is assumed to 
be 1000. Using this value for specific power consumption 
of watts/watt cooling capacity, the required capacity in 
watts for the helium refrigerator can be determined. ^ 

Capacity Weight Volume 

10 watts 2000 lb. 100 ft 3 

A refrigeration unit of the above dimensions would 
be used to supply the cooling requirements for one super- 
conducting machine. Four generators and two motors would 
then require the use of six refrigerators of 10 watt 
capacity. Cross-connect capabilities could also be 
included as a backup for each system. If the liquefier/ 
refrigerator system for one of the superconducting machines 
should suffer a breakdown, the appropriate cross- connects 
could be made and the unit in question would be cooled 
by liquid helium from another machine's liquefier/refrigera- 
tor. 
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4.4.1 Cable System 



Marine cable systems used for connecting main 
generators to the main propulsion motors are sized for the 
allowable losses caused by joule effects. Cable weight 
can be a significant component in electrical propulsion 
systems. The cable weight depends on the voltage, current 
and allowable resistance per length. The advantage of 
ac propulsion systems over dc propulsion systems is that 
significantly higher voltages can be used. Higher voltage 
systems require a lighter cable, which is important in 
systems requiring minimum weight and volume propulsion 
machinery. 

Cables for ac propulsion systems have been studied 
to show the relationship between weights and losses. 

These studies indicate that a light-weight cable system 
is possible with air-cooled cables, see Table 4.4. 

The results in Table 4.4 are based on an allowable voltage 
drop of 9*8 V/lOO ft. The area shown in Table 4.4 is 
the total copper cross-section area necessary to carry 
the 3-phase current at the given voltage. A 6900 volt 
system was used to produce the lightest weight cabling 
system. At 6900 volts, the maximum current will be 
1,260 amps and the maximum transmission loss will be 28 KW 
per 100 ft., resulting in a .1 percent power loss. 
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Table 4,4 



Voltage 


Current 


Area 


Weight 


Loss 


4160 V ac 


2.100 A 


1.56 in 2 


6 lb/ft 


45 KW/100 ft 


6900 V ac 


1.260 A 


1.04 in 2 


4 lb/ft 


28 KW/100 ft 



For a 3-pbase ac system, three conductors of the 
combined size and weight given above must be used to 
transmit power from the generators to the motors. This 
works out for a single conductor to be .384 inches in 
diameter and 1.33 lbs/ft. The weight of the conductor 
is further increased by about 40 percent with the addition 
of armor, lead shield and insulation. Each conductor 
is also surrounded by a 4 inch conduit which provides 
forced-air cooling and an added margin of insulation. 

If this conduit were constructed out of aluminum, it 
would weigh approximately .62 lb/ft for each conductor. 
Adding all of the weight together results in a total 
weight for the three conductors of 7*46 lb/ft. 

4.4.2 Switch Gear and Control Units 

There is very little information on full scale 
switch gear such as would be used on shipboard with 
superconducting electric machinery. Most of the techno- 
logical work performed to date has been on laboratory 
size machinery, the results of which can be used to 
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predict the size and weight of the full scale machinery. 

The product of one such study by Rains^ 1 ^ indicates 
that the switch gear will weigh approximately 2.9 tons 
and occupy approximately 600 ft^. The switch gear and 
control units consist of the actual switching mechanism 
which controls the direction of the motor and an electric 
converter which controls the speed of the motor. The 
electronic converter will most likely be of the cyclo- 
converter type. A cycloconverter is a solid state device 
that converts three-phase fixed frequency power to three- 
phase variable frequency power. The frequency of the 
output power from the cycloconverter directly controls 
the speed of the motor. A cycloconverter of sufficient 

power rating to accomodate a 40,000 HP motor will weigh 

3 

approximately 3*9 tons and require 200 ft-'. Two cyclo- 

3 

converters and the switch gear will occupy 1000 ft-' and 
weigh approximately 10.7 tons. This study also indicates 
that there will be 19.2 tons of miscellaneous equipment 
which will be added to cover items omitted from this thesis. 

a 
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CHAPTER 5 

COMPARISON OF ELECTRIC DRIVE WITH 
MECHANICAL DRIVE FOR SHIP PROPULSION 

5.1 Introduction 

The size and weight of the machinery components 
required to convert the propulsion plant to electric 
machinery were calculated in Chapter Four. In this 
chapter, these components will be combined with the re- 
mainder of the propulsion machinery to predict the result- 
ant size and weight of the new electric propulsion system. 
The new system will be compared to the old mechanical 
drive system to measure the overall savings produced by 
the superconducting electric propulsion system. The 
comparison will be made for two different baseline ships 
as described in Chapter One. The first being the Model 
Baseline and the second the DD963 Baseline. 

The Model Baseline is the DD963 as synthesized on 

MM 

1 the ship synthesis model 7 without the excess volume 
present in the current DD963 design. The DD963 Baseline 
is the DD963 as synthesized with all of the excess volume 

I included in the design. A brief description of the Ship 

I ( 14 ) 

i, Synthesis Model for Naval Surface Ships; 7 along with 

I a sample input and program output is shown in Appendix H. 

I Table 5*1 contains the computer synthesized characteristics 
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Table 5.1 



CHARACTERISTICS OF MODEL BASELINE 
AND DD963 BASELINE 



CHARACTERISTIC 

Max. Sustained Speed, knots 
Sustained Speed SHP, horsepower 
Endurance Speed , knots 
Endurance Speed SHP .horsepower 
Length between perpendiculars, ft. (LBP) 
Beam, ft. (B) 

Draft, ft. (T) 

c p 

C x 

VCG Full Load, ft. 

GM/B 

Average Depth, ft. 

Accomodations 
KW Installed 

Full Load Displacement ,tons( a) 

Light Ship Displacement .tons 

Variable Loads Weight, tons 

WTGP1 , tons 

WTGP2 , tons 

WTGP3, tons 

WTGP4 , tons 

WTGP5, tons 

WTGP6 , tons 

WTGP7, tons 

Weight Margin, tons 

Total Internal Volume, cu.ft.(v^) 

Hull Volume.cu.ft. (v^) 

Superstructure Volume, cu. f t. ( v ss ) 

Full Load Ship Density, lbs/cu.ft. 
Military Mission Volume , cu. ft . 
Personnel Volume.cu.ft. 

Ship Ops Volume.cu.ft. 

Payload Volume Fraction (VOL PAY/ v ) 
Personnel Volume Fraction(VOL PERS/ v ) 
Ships Ops Volume Fraction(VOL OPS/ v ) 
Payload Weight Fraction (W PAY/ a ) 
Personnel Weight Fraction (W PERS/ a ) 
Ship Ops Weight Fraction (W OPS/ a ) 



MODEL 


DD 963 


BASELINE 


BASELINE 


33.9 


32.9 


80,000 


80,000 


20 


20 


10,548 


11,483 


529 


529 


54.6 


55.8 


16.8 


18.8 


• 59 


.59 


.83 


.83 


21.96 


22.02 


.10 


.10 


40.03 


40.67 


298 


298 


6000 


6000 


6906 


7885 


4936.7 


5827 


1868.8 


1959 


2465.2 


3137.1 


789.2 


789.2 


282.7 


296.8 


207.4 


250.3 


569.1 


739.8 


499.8 


454.3 


159.2 


159.2 


100 


100 


945,470 


1 , 013,880 


761,020 


772,581 


184,449 


241,298 


16.36 


17.42 


110,753 


159,298 


251,990 


251,990 


582,728 


602,592 


.12 


.16 


. -27 


.25 


.62 


.59 


.05 


.05 


.04 


.03 


.47 


.44 
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for each of these two ships. 

In Table 5.1, Weight Groups 2 and 7 (WTGP2 and WTGP7) 
are the same for both ships, because the same propulsion 
plant (WTGP2) in type and horsepower and the same armament 
(WTGP7) was specified for both ships. The volume (172,736 
cu.ft.) and weight (789.2 tons) for WTGP2 will be used 
as a baseline for estimating the new volume and weight of 
WTGP2 with superconducting electric machinery. WTGP2 is 
only one of several groups that comprise the entire Machinery 
System. The remaining groups must also be examined for 
changes resululting from a conversion to an electric 
propulsion plant. See Appendix I for a complete listing 
of the contents of each weight group. 

The Machinery System is composed of the following: 
Machinery Box 
Uptakes 

Shafting, Bearings and Propellers 

Maneuvering 

Ventilation 

The Machinery Box can be further broken down into the 
following groups: 

I WEIGHT GROUP 2 -PROPULSION 

Propulsion Units 
Combustion Air Supply 
Propulsion Control Equipment 
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WEIGHT GROUP 2-PROPULSION (cont'd) 

Fuel Oil Service Systems 
Lubricating Oil System 
Propulsion Operating Fluids 

WEIGHT GROUP 3-ELECTRIC PLANT 

Electric Power Generation 
Power Distribution Switchboards 
Electric Power Generator Fluids 

WEIGHT GROUP 5-AUXILIARY SYSTEMS 

Air-Conditioning Systems 
Refrigerating Spaces, Plant & Equipment 
Aviation Fuel & Lube Oil System, Sewage System 
Compressed Air System 

Auxiliary Steam, Exhaust Steam & Steam Drains 
Distilling Plant 

Auxiliary System Operating Fluids 
WEIGHT GROUP 6-OUTFIT AND FURNISHINGS 
Ladders & Gratings 

Weight Groups 3,5 and 6 will not be directly affected by 
the introduction of an electric propulsion system as none 
of the machinery components for these groups will be 
physically changed or replaced. The only changes will 
come from secondary effects as a result of the physical 
dimensions of the ship being modified to accomodate a 
smaller and lighter weight propulsion system. 

The uptakes will change is weight and volume only if 
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the number or size of propulsion unit gas turbines are 
changed, or if electric propulsion allows a major rearrange 
ment of the machinery spaces. The electric plant gas 
turbine's requirements are also considered when sizing 
the uptakes. 

The propellers, shafting and bearings will be 
directly affected by the modification to the electric 
propulsion. The Controllable-Reversable Pitch (CRP) 
propeller will be replaced by a lighter fixed pitch pro- 
peller. Shafting and bearings will be greatly reduced, 
as the electric motors will be coupled to the propellers 
through significantly shorter shafts. The majority of 
the weight for shafting and propellers and all of the 
volume for separate shaft alleys will be eliminated. 

Maneuvering is comprised of steering systems and 
rudders which will not change for either ship when modi- 
fied to electric propulsion. Maneuvering is a function 
of ship size and will vary only slightly as ship volume 
and displacement are changed when the propulsion system 
is modified. 

Table 5*2 shows the volume and weight breakdown for 
the machinery system. All major changes in the volume 
and weight of the machinery system will come from changes 
in the machinery box, uptakes, shafting, bearings and 
propellers. The majority of the weight and volume changes 
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Table 5.2 



VOLUME AND WEIGHT OF MACHINERY SYSTEM 
FOR MODEL BASELINE AND DD963 BASELINE 



MODEL DD963 

CHARACTERISTIC . BASELINE BASELINE 

Machinery System Volume, cu.ft. 293*597 292,801 

Machinery Box Volume, cu.ft. 195*955 195*955 

Uptakes Volume, cu.ft. 49,150 49,150 

Shafting , Bear ings & Propellers.cu.ft. 3*848 3,848 

Maneuvering Volume, cu.ft. 6,996 6,996 

Ventilation Volume, cu.ft. 37,648 36,852 

VOL MACH SYS A *31 .29 

VOL MACH BOX/V .21 .19 

VOL UPTAKES/v *05 .05 

VOL Shafting, Bear ings&Propellers/ v .004 .004 

VOL Maneuver ing/v .007 .007 

VOL Ventilation/ v .036 .036 

Machinery System Weight, tons 1259.1 1205.7 

Machinery Box Weight, tons 722.2 697*3 

Uptakes Weight, tons 130.5 130.5 

Shafting, Bearing&Propellers Weight , tons253*l 253*1 

Maneuvering Weight, tons 81.7 66.0 

Ventilation Weight, tons 71.6 58.7 

Weight Mach Sys/DISPLACEMENT .16 .18 

Weight Mach Box/DISPLACEMENT .090 .10 

Weight Uptakes/ DISPLACEMENT .017 .019 

Weight Shafting, Bear ings&Propellers/ 

DISPLACEMENT .033 .037 

Weight Maneuvering/DISPLACEMENT .011 .010 

Weight Ventilation/DISPLACEMENT .009 .009 
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in the machinery box will result from WTGP2 variations. 

5.2 Conversion to an Electric Propulsion Machinery System 

The results of Chapter 4 sire used to size the machinery 
system by modifying the propulsion system with superconduct- 
ing electric machines. The new weight and volume of the 

machinery system were used as input to the ship synthesis 
(14) 

model to determine the impact on the overall ship 

system. The preliminary output - from the synthesis pro- 
cess indicated that a baseline ship with four gas turbines 
converted directly to electric drive would not meet the 
no-change- in- payload, sustained speed or endurance range 
requirements as set forth in Chapter 1. 

When the volume for the propulsion system went down, 
the volume for the payload went up because the size of 
the ship remained constant. Extra volume for the payload 
is unuseable when a ship is weight limited. More payload 
in the form of extra weight cannot really be added to take 
advantage of this extra volume. If this volume were to 
remain in the ship, it would only serve to make the pre- 
sent configuration of the ship space inefficient and waste- 
ful. 

The sustained speed could be met within the tolerances 
of the ship synthesis model itself. The top speed of the 
modified ships varied by less than 2 %. A 2$ change can 
be considered to essentially meet the no-increase in top 



speed criteria. 

By far, the largest deviation from the baseline 
characteristics was in the endurance range, which increased 
up to 25$. (Section 4.3 contains the calculations for the 
endurance range.) A range increase of this amount was 
unsatisfactory when compared to the inital criteria of 
no-range change allowed in the modified ship. To keep 
the range the same, over 200 tons of fuel can be removed. 

To meet the no-change in mission performance restrictions 
on the modified ship, a smaller three-engined ship with 
less horsepower was then investigated to determine if it 
could take advantage of this fuel weight savings. 

The results from Chapter 4 are summed up in Table 5*3 
for the four-gas turbine, four- generator , two-motor 
propulsion system and the three-gas turbine, three-generator, 
two-motor propulsion system. The four- engined propulsion 
system at 504.1 tons is 36 $ lighter than the baseline 
propulsion system of 789.2 tons. A 49$ weight savings 
is realized by the 401.2 ton three-engined propulsion 
system. The volume occupied by Group 2 propulsion 
machinery is not given in any references. Reference 10 
does contain machinery room arrangement plans without 
specific dimensions being listed. The specific dimensions 
of the gas turbines are listed and could be used as a 
reference for measuring the dimension^ of the machinery 
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Table 5.3 

WTGP2 WEIGHTS FOR BASELINE SHIP. 
FOUR -ENG I NED AND THREE-ENGINED ELECTRIC 
PROPULSION SYSTEMS 



CHARACTERISTIC 



4-ENG. 3-ENG. 
BASELINE SHIP SHIP 



Gas Turbines 




81.25 


60.9 


Generators 




16.0 


12.0 


Motors 




30.2 


23 


Cryogenic System 




5.36 


4.5 


Foundations 




12.3 


9.3 


PROPULSION UNITS. TONS 


244.14 


145.11 


109.7 


Shafting & Bearings 




13.6 


13.6 


Propellers 




36 


36 


Cabling (400 ft.) 




1.33 


1.0 


SHAFTING , BEARI NGS&PROPELLERS , 


TONS 253.1 


5 Q...93 


i>o.6 


Combustion Air System, tons 


58-3 




43.8 


Uptakes, tons 


130,5 


130,5 


98.1 


Switch Gear 




10.7 


8.5 


PROPULSION CONTROL, TONS 


io -.97 


19.2 


15.55 


Fuel Oil Service System, tons 


10.1 


10.1 


ZA 


Lube Oil System, tons 


31.2 


20.0 


16.1 


Repair Parts, tons 


8^5 


JL-5 




Operating Fluids, tons 


42.2 


42.2 


37.1 


Equipment (miscellaneous) 




19.2 


14.4 


TOTAL WTGP2 


789.2 


504.1 


401.2 


WTGP2 VOLUME CU. FT. 


172,736 


143 , 946 


115.157 


WTGP2 VOLUME REDUCTION £ 




16 


33 
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spaces and the equipments contained within them. 

By this method of measure, Group 2 volume was 
determined to be 172,736 cu.ft. The volume of the mechani- 
cal drive machinery removed, such as reduction gears, was 
subtracted from this value. The volume required for the 
superconducting electrical drive machinery was added to 
obtain the values of 143,946 cu.ft. and 115,157 cu.ft. 
for the four-engined and three-engined ships respectively, 
resulting in a propulsion plant space reduction of 1 (Sfo 
for the four- engined ship and 33$ for "the three- engined 
ship. More space could be saved by a complete rearrange- 
ment of the machinery in each engine room, as the location 
of the propulsion gas turbines is no longer limited to 
one specific location by the restrictions of the reduction 
gear and propeller shaft. The result of a simple substi- 
tution of the superconducting electric generators for the 
reduction gears and shafting is shown in Figs. 5»1 and 
5.2. The shortening of each engineroom by 8 feet, as 
shown in the machinery arrangement drawings, plus the 
addition of the two electric motor rooms(12x8xl0 ft.) 
produces the 16^ space reduction given in Table 5.3* 

The value of 115,157 cu.ft. for the three-engined ship 
was obtained in the same manner. A sample machinery 
space arrangement for the DD 963 baseline and the DD 963 
electric drive ships is shown in Fig. 5*3» 



\ 
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ENGINE ROOM NO. 2 SIZE AND MACHINERY ARRANGEMENT FOR 
BASELINE SHIP AND MODIFIED WITH SUPERCONDUCTING ELECTRIC MACHINERY 



DD 963 BASELINE 



Machinery Space Volume occupies 30$ of total 
enclosed volume. 



- Machinery and Uptakes 
g§g§ - Fuel Storage 




Machinery Space Volume occupies 25$ of total 
enclosed volume. 



MACHINERY SPACE ARRANGEMENTS FOR BASELINE SHIP 
AND ONE OF THE POSSIBLE ARRANGEMENTS FOR A SHIP 
MODIFIED TO ELECTRIC PROPULSION 



Figu r e ^ .3 
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5.3 Comparison of Baseline and Modified Ships 



Machinery System weights and volumes derived in 
Section 5*2 for the electric propulsion system were used 
as input into the ship synthesis model. The calculated 
modified ship characteristics were compared to the base- 
line characteristics. The characteristics of the Model 
Baseline, Model Electric propulsion with four-engines 
and the Model Electric propulsion with three- engines are 
given in Table with the same characteristics for 

the DD963 in Table 5*6. Tables 5*5 and 5*7 contain the 
weight and volume breakdown for the Model and DD963 ships 
based on a functional grouping of ship functions. 

5.3.1 Analysis of Model Ships 

When the baseline and the four-engine electric 
drive are compared, there are three major points of 
variance: the range, WTGP2, and the weight of WTGP2 in 
pounds divided by the shaft horsepower. The large differ- 
ence in range can be attributed to the manner in which the 
two different plants are operated. For cruising at the 
endurance speed of 20 knots, the mechanical drive ship 
is required to operate two gas turbines, one for each 
shaft. The endurance shaft horsepower is 10,5^8 HP which 
works out to approximately 5300 HP per gas turbine. Fig. 2. 3 
in Chapter 2 indicates that the specific fuel consumption 
(SFC) of two engines, each operating at 5300 HP, is 
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Table 5.4 



CHARACTERISTICS OF MODEL BASELINE AND MODEL ELECTRIC 
PROPULSION SYSTEM SHIPS 



LBF 


MODEL 

BASELINE 

529 


B 


54.6 


T 


16.8 


SHF INSTALLED 


80000 


SHF END 


10548 


V fp 


33-9 


V end 


20.0 


Range 


6000 


Displacement 


6906 


WTGP1 


2465.2 


WTGP 2 


789.2 


WTGP3 


282.7 


WTC-P4 


207.4 


WTGP5 


569.1 


WTGP 6 


499.8 


WTGP7 


159.2 


Loads 


I 869.8 


l'JT Margin 


100.0 


v Total 


945470 


WT PAY/a 


.05 


WT PERS/a 


.04 


WT OPS/ A 


.47 


VOL PAY/v 


.12 


VOL PERS/v 


.27 


VOL OPS/ v 


.62 


WTGP2/SHP 


22.1 


VOL MACH BOX/SHP 


2.45 


VOL HAB/MAN 


724 


SHP/a 


11.58 


VOL MACH SYS/V 


.31 


WTGP2/A 


.11 


VOL MACH BOX 


195955 


VOL UPTAKES 


49150 


VCL SHAFT. &BEAR. 


3848 



MODEL ELECTRIC MODEL ELECTRIC 
4 ENGINES 3 ENGINES 



529 


501.3 


54.6 


53.7 


16.8 


16.1 


80000 


60000 


10332 


9491 


34.3 


33.4 


20.0 


20.0 


7874.3 


6000 


6657.4 


5938.3 


2465.2 


2181.8 


504.1 


401.2 


282.7 


265.2 


207.4 


207.4 


569.1 


561.3 


499.8 


480.8 


159.2 


159.2 


1869.8 


1771.2 


100.0 


100.0 


945470 


869832 


.05 


.064 


.04 


.04 


.44 


.45 


.15 


.13 


.27 


.30 


.59 


.58 


14.1 


15.0 


2.16 


2.49 


724 


724 


12.0 


10.1 


.28 


.26 


.076 


.068 


172675 


149253 


48789 


34973 


0 


0 
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Table 5.5 

FUNCTIONAL GRO UP W EIGHTS AND VOLUMES OF MODEL 
BASELINE AND ELECTRIC PROPULSION SYSTEM 





MODEL 

BASELINE 


MODEL ELECTRIC 
4 ENGINES 


MODEL ELETRIC 
3 ENGINES 


Vol. Hull 


761020 


761020 


683734 


Vol. Superstructure 


184449 


184449 


168802 


Vol. Total 


945470 


945470 


852536 


VOLUME 


Military Mission 


110753 


139038 


110753 


Personnel 


251990 


251990 


251990 


C ontrol 


56758 


56758 


52975 


Mach. Sys. 


293597 


265312 


221313 


Deck Aux. 


5194 


5194 


4662 


Maintenance 


18614 


18614 


17727 


Stowage 


69736 


69736 


58918 


Tankage 


21225 


21225 


19466 


Pass & Access 


117605 


117605 


114732 


WEIGHT 


Military Mission 


379.8 


379.8 


379.8 


Personnel 


247.8 


247.8 


247.8 


Control 


81. 3 


81. 3 


78.1 


Mach. Sys. 


1205.7 


930.6 


818. 5 


Deck Aux. 


97.9 


97.9 


97.9 


Maintenance 


101.9 


101.9 


100.1 


Stowage 


1732.4 


1732.4 


1277.4 


Pass & Access 


14.9 


14.9 


14.2 


Hull Group 


2502.3 


2502.3 


2181.8 


Ship Sys. 


551.4 


551.4 


529.3 
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Table 5.6 

CHARACTERISTICS OF DD963 BASELINE AND DD963 ELECTRIC 



PROPULSION SYSTEM SHIPS 



DD 963 

BASELINE 



LBP 


529 


B 


55.8 


T 


18.8 


SHP INSTALLED 


80000 


SHP ENDURANCE 


11483 


V sus 


32.9 


V end 


20 


Range 


6000 


Displaement 


7885.0 


WTGP1 


3137.1 


WTGP2 


789.2 


WTGP3 


296.8 


WTGP4 


250.3 


WTGP5 


739.8 


WTGP 6 


454.3 


WTGP7 


159.2 


Loads 


1958.6 


Wt. Margin 


100 


v Total 


1013880 


WT PAY/a 


.05 


WT PERS/a 


.03 


WT OPS/a 


.44 


VOL PAY/v 


.16 


VOL PERS/v 


.25 


VOL OPS/v 


• 59 


WTGP 2 /SHP 


22.1 


VOL MACH BOX/SHP 


2.45 


VOL HAB/MAN 


724 


SHP/a 


10.14 


VOL MACH SYS/v 


.29 


WTGP 2 /a 


.10 


VOL MACH BOX 


195955 


VOL UPTAKES 


49150 


VOL SHAFT & BEARINGS 


3848 



DD 963 ELECTRIC 
4 ENGINES 


DD 963 ELECTRIC 
3 ENGINES 


529 


509.8 


55-8 


54.7 


18.8 


17.18 


80000 


60000 


11237 


10418 


33-2 


32.26 


20 


20 


7848 . 1 


6000 


7568.4 


68 00.7 


3105.6 


2757.7 


504.1 


401.2 


296.8 


275.9 


250.3 


250.3 


739.8 


735-5 


454.3 


445.9 


159.2 


159.2 


1958.6 


1864.6 


100 


100 


1013880 


912864.9 


.05 


.06 


•03 


.035 


.41 


.43 


.18 


.17 


.25 


.28 


• 57 


• 55 


14.1 


15.0 


2.16 


2.49 


724 


724 


IO .52 


8.82 


.26 


.24 


.07 


.060 


172675 


149253 


48789 


34973 


0 


0 
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Table 5.7 



FUNCTIONAL GROUP WEIGHTS AND VOLUMES OF DD963 
BASELINE AND ELECTRIC PROPULSION SYSTEM 





DD963 

BASELINE 


DD963 ELECTRIC 
4 ENGINES 


DD963 ELECTRIC 
3 ENGINES 


Vol. Hull 


772581 


772581 


707468 


Vol. Superstructure 


241298 


241298 


205396 


Vol, Total 


1013880 


1013880 


912865 


VOLUME 


Military Mission 


159298 


186787 


159298 


Personnel 


251990 


251990 


251990 


Control 


60633 


60633 


56912 


Mach. Sys. 


292801 


265312 


221013 


Deck Aux. 


5194 


5194 


4817 


Maintenance 


21330 

72459 


21330 

72459 


18 003 


Stowage 


64573 


Tankage 


28171 


28 171 


20808 


Pass. A Access 


122005 


122005 


115451 


WEIGHT 


Military Mission 


392.8 


392.8 


392.8 


Personnel 


238.8 


238.8 


238.8 


Control 


120.0 


120.0 


119.4 


Mach. Sys. 


1259.1 


974.0 


856.3 


Deck Aux. 


115.1 


115.1 


115.1 


Maintenance 


92.2 


92.2 


90.7 


Stowage 


1806.1 


1806.1 


1427.4 


Tankage 


0 


0 


0 


Pass. & Access 


13.5 


13.5 


12.8 


Hull Group 


3126.1 


3126.1 


2790.4 


Ship Sys. 


602.3 


602.3 


578.2 
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.64 lbs/hp-hr. The synthesis model calculated that for 
a mechanical drive ship with an SFC of .64 requires 1606 
tons of fuel for an endurance range of 6000 nautical miles. 
At 6784 lbs/hr fuel consumed for 300 hours (time required 
to cover 6000 mi. at 20 knots) the ship would consume 
908 tons of fuel with the remaining 698 tons of fuel 
being consumed to provide heat and electricity. Cruising 
at endurance speed requires only one gas turbine operating 
in an electric drive ship. Fig. 2.5 in Chapter 2 shows 
the SFC of one engine operating at 10,548 hp to be 
.5 lbs/hp-hr. When the specific fuel consumption drops 
to .5 lbs/hp-hr with an endurance SHP of 10,548 hp, it 
requires 393 hours to consume 908 tons of fuel. 

Cruising for 393 hours at 20 knots produces an 
endurance range of 787^ nautical miles, which is consider- 
ably greater than the maximum required range of 6000 miles. 
The endurance range can be reduced to the desired 6000 miles 
by removing 202 tons of fuel for a 13$ fuel reduction. 

This reduces the full load displacement to 6453 tons. 

For this displacement the ship is unstable due to the 
decreased draft. More weight would have to be added to 
make the ship stable again eliminating any advantages 
received from reducing the fuel weight. To produce a 
stable ship requires a reduction in the size of the hull, 
reducing the buoyancy and increasing the draft. A lighter 
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and smaller ship needs less power to provide the required 
full speed of 30+ knots; therefore, the primary reason 
for the three-engined propulsion plant ship. A smaller 
ship automatically implies less fuel (greatly reducing 
fuel weight) for the same endurance range. The smaller 
ship with an endurance horsepower of 9. **91 and an SFC of 
.5 requires only 636 tons of propulsion fuel to cover 
6000 nautical miles. This is approximately a JQrfo savings 
in propulsion fuel (a 1?^ fuel saving overall) for a ship 
that can perform identical functions at the same efficiency 
as the larger ship. The percent change from the baseline 
ship to the electric drive ships are shown in Table 5 * 8 . 
Only those items which vary from the baseline values are 
listed. For example, the weight of armament does not 
change nor does the weight for Military Mission; therefore, 
they need not be listed in Table 5 * 8 . 

An examination of Table 5*8 indicates why the three- 
engined ship is preferred over the four- engined ship. 

The four-engined ship has large increases (20 ft or more) 
in range, volume, military mission and the payload volume 
fraction resulting in wasted space. The three-engined 
ship had a 20fo increase only in the payload weight fraction 
indicating a more efficient ship design. The primary 
function of a Navy ship is to deliver as much military 
payload as possible where it is needed. Provided the 
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Table 5.8 

CHANGE FROM MODEL BASELINE 
PERCENT CHANGE FROM THE BASELINE TO ELECTRIC DRIVE, 
LISTING ONLY THOSE ITEMS WHICH CHANGE SIGNIGICANTLY 
FROM THE MODEL BASELINE SHIP 



4-ENGINE ELECTRIC 3-ENGINE ELECTRIC 
% CHANGE* FROM ?£ CHANGE* FROM 

BASELIHE BASELINE 



LBP 

SHP INSTALLED 

V fp 

Range 

Displacement 

WTGP1 

WTGP2 

Volume Total 

Volume Military Mission 

Volume Control 

Volume Machinery System 

Volume Stowage (Includes Fuel) 

Volume Tankage 

Weight Machinery System 

Weight Stowage (Includes Fuel) 

Weight Hull Group 

Weight Payload/Displacement 
Weight Personnel/Displacement 
Weight Operations/Displacement 
Volume Payload/Volume Total 
Volume Personnel/Volume Total 
Volume Operations/Volume Total 
WTGP2/Shaft Horsepower 
Volume Mach. Box/Shaft Horsepower 
Shaft Horsepower/Displacement 
Volume Mach. System/Volume Total 
WTGP2/t)isplacement 

Volume Machinery Box 

Volume Uptakes 

Volume Shafting & Bearings 

Full Load Ship Density lbs/cu.ft. 



0 

0 

- 1.2 

-31.2 

3.6 
0 

36.1 

0 

-25.5 

0 

9.63 

0 

0 

22.8 
0 

0 

0 

0 

6.4 

-25.0 

0 

4.8 
36.2 

11.8 
-3.6 

9.7 

30.9 

11.9 
.7 
100 

3.6 



5.2 
25 

1.5 

0 

14.0 

11.5 

49.1 
8 

0 

6.7 

24.6 

15.5 

8.3 

32.1 

26.3 
12.8 

- 20.0 

0 

4.2 

-8.3 

- 11.1 

6.5 
32.1 
- 1.6 
12.8 
16.1 

38.2 

23.9 
28 • 8 
100 

3.5 



* (-) indicates an increase from the baseline 
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payload is the same and the ships have the same speed 
and endurance characteristics, a ship that has a payload 
which is 6$ of its total weight is more efficient than 
one which has a payload of 5$ of its total weight. 

The overall summation of Table 5*8 is that a three- 
engined electric drive ship which is 14$ lighter, 8$ 
smaller, has a 24$ smaller and 50 $ lighter propulsion 
system and carries 17$ less fuel, can deliver an identical 
payload at the same speed and range as the larger and 
more expensive mechanical drive baseline ship. 

5.3.2 Analysis of DD963 Ships 

As for the four-engined model electric drive, the 
DD 963 electric drive four-engined ship has a much greater 
range than the baseline ship. For this reason, the three- 
engined ship was synthesized to eliminate the excesses 
found in the four- engined ship. The DD 963 three- engined 
electric drive ship also has a 20$ increase in the payload 
weight fraction indicating a more efficient ship design. 

The overall summary of Table 5*9 for the DD 963 is a 
three-engined electric drive ship which is 14$ lighter, 

10$ smaller, has a 25 $ smaller and 50 $ lighter propulsion 
system and carries 17 $ less fuel, can carry an identical 
payload at the same speed and range as the larger baseline 
ship. 

5.4 Final Comparison of Model and DD963 
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Table 5.9 

CHANGE FROM DD963 BASELINE 
PERCENT CHANGE FROM THE BASELINE TO ELECTRIC DRIVE, 
LISTING ONLY THOSE ITEMS WHICH CHANGE SIGNIFICANTLY 
FROM THE DD963 BASELINE SHIP 



4-ENGINE ELECTRIC 
# CHANGE* FROM 
BASELINE 



3-ENGINE ELECTRIC 
# CHANGE* FROM 
BASELINE 



LBP 

SHP INSTALLED 

V fp 

Range 

Displacement 

WTGP1 

WTGP2 

Volume Total 

Volume Military Mission 

Volume Control 

Volume Machinery System 

Volume Stowage (Includes Fuel) 

Volume Tankage 

Weight Machinery System 

Weight Stowage (Includes Fuel) 

Weight Hull Group 

Weight Payload/Displacement 
Weight Personnel/Displacement 
Weight 0 per at i ons/D i spl ac e ment 
Volume Payload/Volume Total 
Volume Personnel/Volume Total 
Volume Operations/Volume Total 
WTGP2/Shaft Horsepower 
Volume Mach. Box/Shaft Horsepower 
Shaft Horsepower/Displacement 
Volume Mach. System/Volume Total 
WTGP2/Displacement 

Volume Machinery Box 

Volume Uptakes 

Volume Shafting & Bearings 

Full Load Ship Density 



0 


3.6 


0 


25 


-.9 


1.9 


30.8 


0 


4.0 


13.8 


1 


11.2 


36.1 


49.1 


0 


9.9 


-17.3 


0 


0 


6.1 


9.63 


24.6 


0 


10.8 


0 


5-4 


22.6 


32.0 


0 


21 


0 


10.74 


0 


-20.0 


0 


0 


6.8 


2.3 


-12.5 


- 6.3 


0 


-12.0 


3.4 


6.8 


36.2 


32.1 


11.8 


-1.6 


- 3.8 


13.0 


10.3 


17.2 


30.0 


40.0 


11.9 


23.8 


.7 


28.8 


100 


100 


4.0 


1.5 



*(-) indicates an increase from the baseline 
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The primary purpose of synthesizing both the super- 
conducting powered Model and DD963 was to use the Model 
as a control case to check on the accuracy of the DD963 
conclusions. The values associated with the model 
synthesis output can reasonably be assumed to be unbiased 
representations of actual changes brought about by the 
conversion to superconducting electric propulsion machinery. 
With the exception of the machinery box weight and volume, 
the Model weights and volumes for all ship functions are 
generated within accepted design lanes based on past design 
practices and philosophies. In the Model, there is no 
wasted or excess weight and volume which can be removed 
to provide a false indication of realized space and weight 
savings when conversion to electric propulsion is made. 

The DD963 is designed outside some of the accepted design 
lanes and consequently has excess volume which must be 
accounted for when analyzing the savings brought about 
by the conversion to superconducting electric machinery. 

The variation of each of these ships from their respective 
baseline ships is shown in Table 5«10» Those characteris- 
tics which are in disagreement are listed in the upper 
half of the table while the characteristics which agree 
are listed in the lower half. 

The small percentage change in the full load ship 
density for the DD963 indicates that much of the excess 
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Table 5.10 



CORRELATION OF MODEL AND DD963 
PERCENT CHANGE OF 3-ENGINED MODEL AND 3-ENGINED DD 963 
FROM THEIR RESPECTIVE BASELINE SHIPS, SIGNIFICANT 
DIFFERENCES ONLY ARE LISTED 

MODEL DD963 

% CHANGE*FROM % CHANGE* FROM 

BASELINE BASELINE 

CHARACTERISTICS WITH VARIANCE 



LBP 


5.2 


3.6 


V fp 


1.5 


1.9 


Volume Total 


8 


9.9 


Volume Stowage 


15.5 


10.8 


Volume Tankage 


8.3 


5.4 


Weight Stowage 


26.3 


21 


Weight Hull Group 


12.8 


10.74 


Weight Operations/Displacement 


4.2 


2.3 


Volume Payload/Volume Total 


-8.3 


-6.3 


Volume Mach. System/Volume Total 


16.1 


17.2 


WTGP2/fc)isplacement 


38.2 


40.0 


Full Load Ship Density 


3.5 


1.5 


CHARACTERISTICS WITHOUT VARIANCE 




Displacement 


14.0 


13.8 


Volume Machinery System 


24.6 


24.6 


Weight Mach. System/Displacement 


32.1 


31 « 9 


Weight Payload/Displacement 


-20.0 


-20.0 


Weight Personnel/Displacement 


0 


0 


Volume Personnel/Total Volume 


-11.1 


-12.0 


WTGP2/Shaft Horsepower 


32.1 


32.1 


Volume Mach. Box/Shaft Horsepower 


-1.6 


-1.6 


Shaft Horsepower/Displacement 


12.8 


13.0 



* (-) indicates an increase from the baseline 
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volume had been removed during the synthesis process, the 
ship would have been much more dense. As it is, the density 
of the Model increased more than the density of the DD963. 

By comparing the changes in total volume, the DD963 only 
lost 2# more of its original baseline volume than did 
the Model. These two items considered together, would 
indicate that very little of the excess volume in the 
DD963 design was lost when the ship was converted to an 
electric propulsion system. 

The smaller decrease in the volumes of stowage and 
tankage for the DD963 follows the conclusion that the 
DD963 is a less dense ship. These two volumes were 
greater to begin with, and when the volume for the fuel 
savings was removed it created a smaller percentage 
change than in the Model ship. Volume is also the cause 
for the full load displacement of the DD963 decreasing 
less than it did for the Model. The displacement of the 
DD963 is forced higher than it need be in order to have 
a hull large enough to enclose the required volume. If 
the volume were allowed to decrease more, the changes in 
full load displacement would be much closer. 

Taking all of the above into consideration, the 
following conclusions can be drawn for powering a ship 
| with superconducting electric machinery: When compared 

to a four-engined mechanical driven ship the three-engined 
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electrical drive ship willa 



be 9 $ smaller in volume 

be 14$ lighter 

carry an identical payload 

have the same endurance and top speed 

have identical habitability standards 

have a 32 $ lighter machinery system 

have a 25 $ smaller machinery system 

use and carry approximately 17 $ less fuel 

have a 50$ lighter propulsion system (WTGP2) 

have a 33$ smaller propulsion system (WTGP2) 
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CHAPTER 6 
Cost Analysis 



6.1 Introduction 

At this time, it would be very difficult to put a 
price on the acquisition and installation of a superconduct* 
ing electric propulsion system. This thesis covers an 
economic comparison of how much such a propulsion system 
could cost and still be economically feasible. The cost 
of the mechanical drive machinery removed will be figured 
plus the difference in operating cost over the life of the 
ship. The operating costs of an electric drive ship are 
based on a 20 year life cycle, due to the fact the DD 963 
is designed on a life cycle of 20 years. 

6.2 Cost of Removed Machinery 

The single largest weight removed from the propulsion 
system is the shafting, bearings and propellers. 191*3 tons 
of shafting and bearings and 10.88 tons of propeller were 
removed. The cost of this removed weight can be assumed 
to be $2000/ton, ^"which results in a $404,360 cost 
savings. 32.4 tons of uptakes at $1000/ton^^ were 
also removed for a cost reduction of $32,400. 

The reduction gears were the last large weight and 
volume pieces of machinery removed. The cost of the reduce 
tion gears is based on the total horsepower rating of the 
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gear and not on the weight of the machinery, as were the 
previous items. A representative cost of the reduction 
gears as installed in the DD 963 can be based on a $20/hp^ 1 ®^ 
figure. At $ 20 /hp, the cost of removing 80,000 shaft 
horsepower of reduction gears is $ 1.6 million. The removal 
of the above mentioned equipments results in $2 million 
savings . 

6.3 Propulsion Plant Operating Costs 

An operating profile of 3<$ underway time per year 
will be chosen as representative of this class of ship. 

While underway, 94# of the time will be at endurance speed 
and 6 # of the time will be at full power. Underway fuel 
consumption will be calculated on a specific fuel consump- 
tion (SFC) of .42 lb/hp-hr while at full power. For the 
electric drive ship SFC at endurance speed is .5 lb/hp-hr, 
while for the mechanical drive endurance SFC is .64 lb/hp-hr. 
The SFC for the ships electrical power generation will be 
.96 lb/kw-hr^ 12 ^for an average 24 hour electric load of 
1600 kw/hr^ 1 ^. The price charged for fuel will bei $16.8/ 
barrel based on 1977 dollars and fuel prices. Total 

( 19 ) 

fuel costs per year are' 



Fuel CostAear =^](SFC i x SHP i x HOURS -AEAR x FUEL COST/lb) 

1 (6.1) 

Manning costs will be considered to be approximately 
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the sane for both the electric drive and the mechanical 
drive ships. For a Navy ship, there should be the same 
number of men of equal pay grades and skills required 
for plant operation and maintenance. The difference in 
life cycle costs for the propulsion plants is calculated 
on acquisition costs and fuel costs. 

6.3.1 Fuel Costs per Year 

(iq) 

Fuel consumed per year is calculated by« v 77 

Fuel/ rnMC ,v = (SFC)( davs underway x 24hr/day x #time at SFC ) 

2240 lb/ton 

(SHP) (7.23 BBLS/TON) (6.2) 

The fuel consumed by the baseline ship is calculated first. 
Fuel(FULL POWER) = ( - 42 ? i 10 9 x ||^ .061(80,000) (7.23) 

= 28,371 bbls 

Fuel (ENDURANCE) = 52,920 BBLS 
Fuel ( ELECTRICAL ) = 12,969 BBLS 
Fuel Consumed/year = 99,670 BBLS 

At $16. 8/BEL, fuel cost for the baseline ship.:*.is 
$1,674,454. 

The fuel consumed by the electric drive ship is* 

Fuel (FULL POWER) = ( »^ 2 ) ( 10 9 * ^ x ° 6 ) ( 60 » ° 9 &) . (7- 2 3) 
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Fuel (Full power)® 12,767 BBLS 
Fuel (Endurance) a 41,344 BBLS 
Fuel (Electrical) = 12,969 BBLS 
Fuel Consumed/year = 67,080 BBLS 



The total annual fuel cost for operating the elctric drive 
ship is $1,126,944. The cost savings in fuel alone is 
$5^7,510 per year. 

6.3.2 P.V. of Life Cycle Plant Costs 

To compute the present value (PV) of plant operating 
costs the discount rate (DR) is assumed to be 6$ (based 
on real value 1977 dollars) . The discount rate factor 
is* 




(1 + DR) L - 1 
DR(1 + DR) L 



(6.3) 



where L is the life of the ship. Present Value (PV) is 



PV = Cost/year (C DR ) 



(6.4) 



The present value of the fuel saved over the 20 years life 
of the ship is* 

C DR = 12 ‘5 

PV = $5^7.510 (12.5) 

-= $6,843,875 

The Present Value of $6,8 million fuel cost savings 
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is valid only if the cost of fuel does not change over 
the next 20 years. One would have to be very naive to 
even assume that the cost of fuel will not continue to 
change. The $6.8 million fuel cost savings is a bottom- 
line figure based on the assumption that fuel costs will 
not rise any faster than will inflation. As the cost of 
fuel continues to rise, the fuel economical ship is even 
more desirable and the fuel savings for the electric 
drive ship continues going up. 

6.4 Limit Cost of Electric Drive 

The total cost of equipment removed plus the present 
value of 20 years fuel savings is 8.48 million. Using this 
value as a guideline, the superconducting electric propulsion 
machinery is economically feasible if the acquisition cost 
of threee generators, two motors, six cryogenic refrigera- 
tors, switch gear and cabling is $8.48 million or less. 

For a 60,000 hp ship, this value breaks down to a cost 
of $141 per horsepower, or $188 per KW. $188 per KW for 
an upper limit cost of a complete superconducting electric 
propulsion system appears to be obtainable utilizing current 
superconducting technology. 

The $8.48 million savings figure does not take into 
account the decreased acquisition cost or the yearly de- 
crease in maintenance costs of a 10$ smaller ship. The 
20 year cost of paint alone will be a substantial savings. 
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Inclusion of these additional savings would only increase 
the upper limit acquisition cost of an electric transmission 
system. 
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CHAPTER 7 



Conclusions and Recommendations 

The results of the design criteria which were applied 
to a propulsion system incorporating superconducting 
electrical machines indicates that development of this 
technology will provide a significant reduction in machinery 
weight and volume. In any comprehensive study of competi- 
tive propulsion systems for a particular ship design, the 
superconducting electric system must be considered a viable 
candidate for increasing the overall ship performance. 

The work associated with this thesis leads to some 
specific conclusions, resulting from the substitution of 
superconducting electric machines for the transmission 
system presently installed in the DD963 and the subsequent 
reduction in ship size and propulsion system size and power. 
Additionally some general conclusions can be drawn for the 
application of these machines in other propulsion systems. 
The specific conclusions can be enumerated: 

1. A 388 LT reduction in machinery plant weight (50$), 
exclusive of fuel, can be realized by the substitution of 
the proposed system for the presently installed. A projec- 
tion of gas turbine fuel consumption improvement allows 
fuel weight to be reduced by 272 LT for a 17# reduction. 
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When machinery and fuel weight savings are considered 
together, there is an absolute weight savings of 660 LT. 

2. Substitution of the proposed system in the DD963 
produced a marked improvement in the volume required for 
the machinery plant, for a volume reduction of 29,000 ft^ 

(1 Qfo) , This reduced volume would still allow sufficient 
space for the performance of machinery maintenance. 

3. A smaller ship carrying an identical payload and 
having identical performance characteristics can be 
constructed to take full advantage of the weight and volume 
savings provided by the proposed propulsion system. The 
construction of the smaller three-engined ship versus the 
original four-engined ship, produced an overall weight 
reduction of 1085 LT for a 14# weight reduction and a 
savings in volume. 

4. A superconducting electric transmission system 
is economically feasible if its acquisition cost is $8.48 
million or less ($141 per horsepower) . 

In addition to the above conclusions, some general, 
observations about the proposed superconducting system 
are: 

1. The use of an electric transmission system allows 
the efficient operation of gas turbines without the problem 
of excessive fuel consumption at off-design conditions. 

The power plant can be divided easily into several 
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generating units, each with a gas turbine and superconduct- 
ing generator, without requiring the individual unite 
being clustered around the reduction gears. Increased 
prime mover dispersion will decrease the vulnerablility 
of the propulsion plant. The number of economical operating 
speeds would correspond to the number of generating units. 

2. Application of the proposed system in any particu- 
lar ship design would not preclude any future changes in 
components that take advantage of technological improve- 
ments. All components of the proposed propulsion system 
are small enough that they could be replaced by an improved 
machine with relative ease. Substitution of original 
equipment with new and improved designs may be very 
desirable over the life of the ship. 

3. The proposed electric drive system has a high 
degree of controllability, primarily due to the ability 
to use all electric controls for power control. 

These conclusions demonstrate the desirability and 
the feasibility of incorporating superconducting electric 
machines in Naval ship propulsion systems. 

This thesis has touched on several areas which require 
5 further investigation before superconducting electric 
'] machines can be used in marine applications. The most 
I important of these is the actual details of machine construc- 
1 tion. In addition to the specific machine design problems, 
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the use of solid state converters for the control of 
synchronous motors and the liquefier/refrigerator oryogenic 
system, as proposed here, will require further exploration 
and development. The poor controllability and torque 
characteristics of synchronous motors at low power levels 
and RPMs would probably necessitate the use of a conventional 
electric motor to start and power the propellers at low 
RPMs. This motor could be in the configuration of a motor 
around the propeller shaft driven by the ships service 
electrical power. Low speed maneuvering in restricted 
waters and emergency propulsion in the event of total 
failure of the superconducting machinery could be provided 
by these conventional motors. 

A detailed engineering design study would have to be 
performed before the proposed propulsion system could be 
substituted for the original in the DD963. Propulsion 
plant machinery rearrangements should be made to prove 
the validity of the conclusions drawn in this thesis. 

Weight and moment calculations would be required to check 
the stability and trim of the modified vessel. 

The preceding suggested areas of development and study 
sure only those which are the most obvious at the conclusion 
of this thesis, and certainly many additional problem areas 
will have to be addressed for theory ±0 be translated into 
reality. 
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APPENDIX A 



Superconducting Machine Design Program 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

14 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

36 

39 

40 

41 

42 

43 

44 

45 

46 



Cl JOB 



SUPERCONDUCTING MACHINE DESIGN ***♦ 








DIMENSION V (8) ,HC(5),JC(5) 










DIHENSION VV (8) ,DV(B) ,VV1 (0) 










DIMENSION CVL(e) 










C 


RFI , TUF, G FK, 


THK , GK A, 


Tll», GAS, 


AJP 




DATA 


V(1), V (2) , V(3), 


V ( 4) , V (5) , 


V (6) , V (7) , 


V(8) 




1/ 


.40, • C26 , .025, 


.05, .02, 


.1 . .02. 


1.28*8 


/ 


DATA 


DV ( 1 ) , DV (2) , DV (3) , 


3V(4), D V (5) , 


DV (6) , DV (7) , 


DV (8) 




1/ 


.04 , .0026, .0025, 


.005, .002, 


.01, .002, 


1.28*7 


/ 


DATA 


D VL ( 1 ) , DVL (2) , 


DVL (3) * DVL (4) 






V 


* 1 , • 1 , 


•1# 


•1 / 






DATA 


DVL (5), DVL (6), 


DVL (7), t VL (8) , TW 






1/ 


• 1 , • 1 , 


• 1 0 


•1, -I / 







D AT % NOHIT 

1 / 15 / 

CFP=CF (V, NOHIT) 

c ***« STEPSIZE DETER HI NITATION FOR OPTI HIZ ATIOM **♦* 

EPSI x 0 .005 
N V =8 
CFO=C f P 

IF(N1I,’1IT . EQ. C) GC TO 100 
DO 2 NN-1, NO HI 7 
DO 5 1=1, NV 
DO 1 11 = 1, NV 

1 VV(II)*V(II) 

T 1 *C P 4V V , NOPIT) 

A X .CC1*VV (I) 

IF (DV (I) » IT * A) DV (I) *VV (I) *. C01 

VV (I)*VV (I) *DV (I) MW 
T2= CF ( VV , N UM IT) 

VV (I)=VV (I ) *DV (I)MW 
T3=CF (VV, NOHIT) 

I MU M3- 2. *12)50,50,53 
50 DX=-CV(I)MW 

IF ( (T1-T2) .IT. C*0) CX = CV (I) *TW 

52 IF( (T1-T2) .EC. 0.0) DX*0 .0 
GO TO 54 

53 DX = DV (I) * (3. *T 1-4.*Y2*Y3) / (2.* Y1-4.* Y2*2.* Y3) *TW ' 

54 DV (I) = V (I) * C VL (I) 

IF (DX .LT. 0.0) DV(I) =V(I)*(-1.0)*DVL(I) 

IF (ABS (CT) .LT. (D VL (I) *V (I) ) ) DV (I)=DX 
5 CONTINUE 

CO 3 KK*1 , N V 

3 VV 1 (KK) -V (UK) ♦TV (KK) 

CF1=CF (VV1 , NOHIT) 

WRITE (6, 75) CF1 , NN 

IF ( (CPO/CF 1) - ( 1 .0*EPSI) ) 56, 56, 55 ; 

55 DO 4 KK=1, MV 

4 V ( KK) =V V 1 ( KK) 

CFO=CF1 

2 CONTINUE 

56 NUMIt-0 i 

CFO»CF(V,NUPIT) 

75 FORH AT ( 1 *,5X,15B COST PONCTION ,5X,E10. 4 , 10X, 14H ITERATION NO., 
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15X.I6) 

47 100 STOP 

48 END 
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<49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 



FUNCTION CP ( V, HUH) 

DIMENSION V (8) ,HC(6) , JC(6) 

DIMENSION H (20) 

REAL LZ , LCA, IA, IAS, LCA , ICK, LCF, LBB, NAAT, NF AT, LAP 

REAL IA, I FI , NA, IN, MSS, HK , MS, MB, KGLOS, JSR, HA 

REAL JC, 12, HP A , KW A , KBL, K8FL, LTH 
REAL KVAPU, NAPA, HU, KK 
REAL KEKL 

OATA HC ( 1) , HC ( 2) , HC(3), RC(4), HC (5) , IIC (6) 

1/ 5 . 5E* 6, 4.38E46, 3.18E+6, 2.0E*6, .796E»6, 0.0 



DATA 


JC{1) 


. JC(2), JC(3), 


JC<4), JC (5) , 


JC (6) 


1/ 


0.0, 


1.0E*8, 2.0E*8, 


3. 0E*fl , 4. 0£* 8, 


4.6B+8 


DATA 


VA, 


PP. PCIE, BPH, AJA 




V 


22.5 B*6, 


1.0, 6.0. 200.0, 3. 5E*6 / 




DATA 


THWAE, 


THWFE, SPA, SPF, NPA, 


KUA, KBL,KBPL, LTH 




1/ 


1.047, 


2.094, . 3, . 5, 3., 


1., 1., .5, .25 


/ 


DATA 


SIGH A A, 


SIGMAK, KVAPtl, VT, 


XT, XI, X 2, PR 




V 


6.0E+7, 


lm 0E+ 7, 1*, 1 . , 


.1, *1# • 1 , .3 


/ 


DATA 


ROSS, 


F, THAX, ROB, 


GKI, GAI, 12, BP 





1/ 8000., 2.0E+11, 4.5E4E, 1800., .02, .02, .05, 1000./ 

DATA ROAL, DM AX , ROPE, GAMMA 

1/ 2600., 2.4E+8, 75CC., 2.4 / 

DATA FCCO, EAL 
1/ 8800., 6.S4E+10 / 

DATA BSHAX, KBKL, PZ 

1/ 1.75, 1.0, 2.65 / 

R P 1 = V (1) 

Til F=V (2) 

GP K= V (3) 

THK^V (4) 

GK ?.= V (5) 

THA =V (6) 

GA S= V (7) 

A J P=V (8) 

PI =3 .14159 



/ 

/ 



MU=PI*4. OE-7 
PSI = A FCOS ( FF) 

P=POl E/2.0 

OMEGA c PPM* P*2. 0* PI/60.0 
HP=V A/746. C*PP 

C**** CA LCT1I ATE MACHINE DIMENSIONS **♦* 

C+*+ * RPO, PKO, RAO A R 10 OTER RADII OF FIELD, DAMPER, AND ARMATURE 
499 R PO = R F I +T11 P 

RKI =F FO*G PK , 

R KO*=R KI*TH K 
RAI-FKO*GKA 
R AO=R A I*T H A 
RS“RAC+GA5 

X=RAI/RAO V 

Y=RrI/RFO 

U=RKC/PS 

W=R AC/3S 



Z=RFC/ES 

ZZ* RKI/R KO 

BA= (FAO+RAI) /2.0 
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93 

04 

95 

96 

97 

98 

99 
100 
101 
102 
103 



104 



105 



106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 

122 

123 

124 

125 
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R F* (RPOtRPl)/2.0 
RSB* (RKC*RKI)/2.0 
C 1 0= T E H A E/ 2 • 0 
Cl 1*THVFE/2.0 
SCI 0* SIN (CIO) 

SC 1 1 = SI N (C 11) 

CMl*CH (P,X,N) 

CS1*CS (P # X * W) 

CC=2 • 0*P 
P P = 2. 0*P 
AA =2 • C-P 

C**** CALC LENGTH OF MACHINE »*** 

C**** Ills INTERNAL BASED REACTANCE **♦* 

XI= (N PA* SORT (2 . )/4. ) * (AJA/AJF)* (SCI 0/SC1 1 ) * ( ( RA O/RPO) ** BB) 
1*(CS1/(CH1*(1.0-Y**BB)))“K«A 
C+++ + PPI IS KVA/nHIT LENGTH BASED ON INTERNAL VOLTAGE 

PPI= (B. / (SOFT (2. ) *FI) ) *OM EG A*M U* AJ A* AJF* (1.0-Y**BE) * (RPO**BB) 
1* (RAO** A A) *CM1 *SC10*SC1 1*KWA*NPA 
C**** LZ IS UNIT LENGTH 
LZ = V A/PPI 

C**** LDA IS EFFECTIVE LENGTH CF ARMATURE 
L D A = (FAO + RAI - ( P FO 4 R FI ) * K BPL) *KEL/P 
PET A=LDA/LZ 

A T- ( (X I ** 2 *B ET A) ♦XI* SIN (PSI) ) / ( 1. 0-XI-*2) 

ALPH A=AT 4 SQFT ( AT **2 + 1 . C 4 2 . C* XI *B ETA *SI N( PSI) ♦ (XI* BETA) **2) 
C**** LA IS ACTIVE ARMATURE LENGTH 
LA= L7* ALPH A 

C**** LAS IS ARMATURE STRAIGHT SECTION LENGTH 
L A S=L A- (RFO 4 RFI) * ( KBFL) * K PL/P 
C**** LOA IS TOTAL EFFECTIVE ARMATURE LENGTH 
LOA = LAM,DA 

C**** TA IS TOTAL P EACTA NCE/U NIT BASED CN INTERNAL VOLTAGE 
X A-X I * (TCA/LA) 

IF (TA . LT. 1.0) GO TO SCO 
TU A = TH A-TH A *0. 1 
GO TC 4S9 

500 V OE= SC RT ( 1 . 0- (XA*COS(PSI) ) **2) -XA*SIN (PSI) 

XD= X A/VO E 

C**** LOK A N C LOF ARE THE EFFECTIVE LENGTH OF CAMPER AND PIELD 
LOK = L A 4 (RK I ♦ PRO) *{KBL/P) 

LO P = L A 4 ( R F I 4 R F C) /P 

C*+** PULS CF THUMP GUESS FOR LENGTR OF EEARING SPAN 
LBR=LCK 4 2. C * IT H 

C**** NA AT IS TOTAL NUHBFR OF ARMATURE AMPERE-TURNS 
N A AT= AJA*THVAE*RAO**2* (1.0-X**2)/2.0 
C**** NFAT IS TOTAL NUMBER OF FIELD AMPERE-TURNS 
NFAT=AJF*THWPE*RFO**2* ( 1 . 0- Y** 2)/2 . 0 
C**** VP T IS GENERATCD V CLT AGE/TU RN/FH ASE 

VPT=32. 0/(2. 0*SQPT(2.0)*PI )*OHEGA*MU*LA*AJF*SC10*SC1 1 
1* (1 . 0-Y **BP) / ( THKAE* ( 1 . C-X *• 2) * R AO** E) * RFO** BD* CM 1 *K »A* VOE 
C 

C**** CALC OF TRANSIENT ELECTRICAL PARAMETERS 
C**** LAF IS ACTIVE ARMATURE LENGTH FOR COUPLING TO CAMPER 
LAK*lAS*(RF0 4 RFI)*KWA* KBKL/P 

C++** KDP AND XDFP APE TRANSIENT AND SUBTHANSIINT REACTANCES 
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TDP* X D* (1. C-4.0*(LA** 2*KW A/ (lCA*LOr) )* (CM (P # X, W) ** 2) / (CS (P,X,W) 
1*CS (F,T,Z) )* (RFO/RAO) **CC* (1. 0-Y* *BB) **2) 

XDPP-TD* (1.C-2.0* ( (LA**2)/(LCA*LCK) ) * (CM 1**2) / (CS (P, X, W) 
1*CS(P # Z7.,U) )* ( ( FKO/RAO) **CC) * ( 1. 0- ZZ*+ B9 ) ** 2) 

C**** TS AND T DPP ARE SHIELD AND SUBTRANSIENT TIME CONSTANTS 
TS*PI*?1U*PKC*(FKC-FKI) ‘SIGHAK/4.0* (1.0*Z7)**CC 
T A IS ARMATURE TIME CONSTANT 

TA = 2. 0* HU^SC 10** 2* SIG M A A* SFA* R AO** 2* CS (P,X,N)/ 

1 ;PI*1HWAF* (1 • 0-X**2) ) *NPA*K«A**2*XDPP/XD 
TD PP-T S* (XDr-XPPF)/(XD-XDPP) 

c 

C**** FIELD enPRENT PISE 

C FIELD CURRENT RISE DURING CRITICAL POST-FAULT SWING 

C XE IS EXTERNAL REACTANCE 

C EFO AND VINF APE INTERNAL AND BUS VOLTAGES 

C IF 1 IS PER UNI 1 n AX FIELD CURRENT 

ia=kvapu/vt 

XE=X1*X2/(X1*X2) +XT 

EFO = SORT (VT**2+ (XD*IA) **2*2. C*VT*XC*I A* SIN (PSI) ) 

VINF=SQRT(VT**2+ (XE*I A) -*2-2.0*VT*XE* I A*SIN (PSI) ) 

DEODARS IN (KVAPU*COS (PSI) *XE/ (VINP*F.FC) ) 

I FI =2. 0* (XD-XDF) *COS (DEO) / (XEP + XE) / (EPO) ♦ 1* 0 
C 

C**** TORQUE TUBE RECUI REHEN I S 

C TORQUE TUBE IS DESIGNED TO CARRY WORST CASE TORQUE FROH 

C A IINE-LINE TEFMINAL FAULT 

C TPP IS PER UNIT WORST CASE TORQUE 

C TT IS WORST CASE TCRQUE 

C SQ IS RADIUS RATIO OF TORQUE TUBE, THIS ROUTINE ATTEMPTS TO FIND 

C THE PROP I R VALUE FCF SQ 

N = 0 

I A*K V APO/V T 

F.PP = SQRT (VT **2 ♦ (XDPP* I A) **2*2. C*VT*X DP P*IA* SI N (PSI) ) 

TPF=1.3*EPP**2/XDPP 

TT~ V A*TPP/CHEGA 

C**** FINDING OPTIMUM INNER FADIUS 

SS= ( (? . 0*P R) /I 6. 0) *ROSS*RFI**2*CHEGA**2 
FR =0* 5 
SQ=0.9 
SO NEW =0 . 0 

C***» WILL TFT 250 TIMES FOR A SOLUTION 

C IF AFTE* THAT MANY TRIES IT HAS NCI FOUND A SOLUTION 

C IT ESTIMATES STRESS FOB AN ALMOST SOLID SHAFT 

1210 IF (N . LT . 25C) GO TO 1220 
51 =2 . 0*SS 

SO=SI* (1,0 -PR) / (3 • 0*PR) 

TADO-2.0*T1/ (PI * RPI ** 3) 

TO= SQFT (TAOO**2*SQ**2) 

TD^SI 

I? (TO .GT. SI) TD=TO 
F PI = C. 0 
GO TC 1300 

C**** REJECT NEGATIVE INNER RADIUS 
C AND INNER FADIUS GREATER THAN OUTER RADIUS 

1220 IF(SQNEW*SQ .LT. O.C) SQNEN*SQ NEW/2 . 0 
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1(0 

1(1 

1(2 

163 

1(4 

1(5 

1(6 

167 

1(8 



1(9 

170 
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172 

173 

174 

175 

176 

177 

178 
17S 
180 
191 
182 
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184 
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186 
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188 
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190 



IF(SCNIW*SQ .GT. .99) SQN I W-SQNFW/2. 0 
SQ«S0*SQNEW 

C**»* SO, SI ARE HALF OF CENTFIFDGAL STRESSES AT OUTER, INNER BA CII 
SO-SS* (2.0* (1.C-PR)/ (3 .0 + FF) *2.0*50* *2) 

SI«SS* (2.0* (1. 0 - FR) / ( 3 . 0 ♦ P F) *SQ**2*2.) 

C**** T AUO AND T MJ I AR E T OR 0 II E STRESSES 

TAUO = 2. 0*TT/(RFI**3*( 1. 0-SC** 4) * PI) 

TAUI = 2.0*SQ*TT/ (RFI**3 * (1 .0-SQ**4 )*PI) 
c **«* T0 AND TI ARE hoHR*S CIRCLE ADDITIONS 
TO^SCRT ( T A 00* * 2 ♦ SO **2 ) 

TI=SQST (TAUI**2*SI**2) 

C**** LARGFF STRESS POINT USED AS CRITERION 
T 1 =T I 

I F (TC .GE. TI) T 1 = TO 

C**** ATTEMPT TO GET WITHIN 95% CF SFECIPIED STRESS, THAT 
IF ( T ?! A X .LT. TI .AND. .95*TWAX .GT. TI) GO TO 1230 
C**** RBI IS SUPFCRI INNER FADIUS 
R0I= SQ* RFI 

C**** TD IS THE STRESS LEVEL USED TO CALCULATE SHAFT STRESS 
C PENALTY FUNCTION 

T D=TK A X 
GO TC 1300 

C**** NEWTCNS METHOD IS USFC TO OBTAIN A NEW GUESS FOB SQ 
C DTI IS RATE CF CHANGE OF STRESS WITH SQ 

C DS AND DT ARE COMPONENT DERIVATIVES 

C SQNEV IS THE ESTIMATE FOB THE CHANGE IN SQ REQUIRED 

1230 IP (TI-TO) 124C, 1259,1250 
1250 N=N*1 

DS=9.0*SQ*SS* ( 1 . 0 ♦PR) / (3.0*PR) 

DTO=T AUO*4. C*SQ**3/ (1 ,0-SC**4) 

DT*TAU0 ♦DTO 

DT1= (SI* DS ♦TAUI*DT)/T1 
SQNEW = FP*(.975*TMAX-T1) /DTI 
GO TC 1210 
1240 N = N* 1 

DS*4 . C * SQ* SS 

DT=T A00*4. C*SQ **3/ (1. 0-SQ**4) 

DTI = (SC*DS+TAUC*DT) /TI 
SQN EW= PR* (,575*TKAX-T1)/DT1 
GO TC 1210 

C**** DAMPEN REQUIREMENTS FOR FAULT CRUSHING LOADS 
C DAM FR IS DESIGNED TO CAFFY WORST CASE CR DS HI NG LOADS 

C FROM A LINE-LINE TERMINAL PAULT 

C STMG IS MAGNETIC STRESS AT SHORT CIRCUIT 

C STB IS TOTAL STRESS FROM A FAULT 

C STCP IS CENTRIFUGAL STRESS AT RATED SPEED 

C DFCF IS CENTRIFUGAL DEFLECTION AT RATED SPEED 

C DFMG IS MAGNETIC DEFLECTIC* FRCK A FAULT 

C DFB IS TOTAL DEFLECTION AT SHORT CIRCUIT 

1300 FT= (RFI*RBl)/2.0 
THSB=RKO-R KI 
WVW-FT/R S 
ZZZ 3 P S F/RS 
C NTH* ABS ( PF) 

SNTH=SQRT ( 1 . 0-PF* PF) ♦ PF/C NTH 
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211 

212 
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214 
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SNDL=XA*CNTH 

CNDL-SQFT (1.0-SNDL*SNDL) 

BA 0 = 4. 2 42* HU /PI* SC 10*B AO* (1.0-X+ (1.0-X**3)/3.0*W**2) *AJA 
ELSS =P I *NU * (1 . €♦ |HSB/RS)**2)/8.0 
NA *NA AT/A J A 

ELA= ( (16.C*LOA*MU* (NA**2) * (SCI 0**2) ) / (P* PI* (SC10**2) 

1* (1.0-X**2)**2)) * SS*K W A** 2 
FLS=El 8*1-5 

EMAP=2.0*m*SC 10* ( PP/RAO) * (I.O-XO. 0/3.0* ( 1 . 0-X**3) * (R AO/RS) * *2) 
1/(TKWAF* (1 .0-X**2) ) 

EMAS=E1AP* (PSB/Rl) 

XDl1=Xt* (1.0-1 .5* Efl AS* EH A S/ELS/ELSS *1 A/LOA ) 

BA 1=E A O/XD 1 1 

B FO =0.66667 *HU /PI* SC 11 *RSB* (RFO/RSB) **3* (1 .0-1 **3) * ( 1. 0-ZZZ** 2) 

1 * A JP 

DELT=ATAN (SNDL /CNDL) 

THET= AT AN (SNTH/CKTB) 

CCC =2. 0*B A 1/(1 .0+7,7Z**2) 

AAA=E AO*S NTH ♦CCC + PFO* CNDL 
BBB=PAO*CNTH+B FO* SNDL 
ROT= SORT (AAA**2*PBB**2*CCC) 

AT3A*A7AN (BBB/AAA) 

EEE=EPO*COS (ATEA-DELT) ♦FAC*SIN (ATBA*THET) 

DDD=EPC*SI N(ATBA-DELT) -BAO*COS (ATE A*TH IT) 

FR 1= (BOT**2-FEE**2-DDD**2)/MU/4.0E4 

PR 2=S0RT ( (FOr**2-EEE**2*DCE**2)**2*(EEE*DDE) **2) 

1*4. O/HU/4 . 0E4 

C**** PR IS MAX RACIAL FORCE AT FAULT 
FR = E f 1 ♦ FP2 

STMB-PR2*2.0*nSB** 2/THSB/TBSB 
STHU=-PR 1*FSB/THSB 
STMG=STHB-STMU 
VLB=RFO**2-PKI **2 
PH SB = VLB*RO AL 

STCF=RRSE*?SE**2*CREGA**2* 1 .CE-4 
STB= SI HB*STCF* STMU 
AIS-RKO**4-RKI**4 
ESB=E A l * AI S 

DFHB=7. 0/3. 0*PF2*2 .0E4/ESE/THSB**3*RSB**4 
DPMU = STrn*FSB/ESE* 1. 0E4 
DPKG= CFMB-CPMU 
DFCP=STCF*Rsn/ESB* 1. 0E4 
DFB=CFMB*DECE*DFMU 
C 

C**** NEGATIVE SEQUENCE LOSSES 
C KK IS SHIELD CURRENT DENSITY 

KK = U.0*AJA-I2*SC10*RKC** (P-1 .0)*RAC** ( 2 . 0-P) * P* BB* CH (P ,X,«)/ 

1 (PI* (1 .0*H**CC) ) 

DDS = SQRT (2. C/ (CHFGA*Hn*SIGHAK) ) 

PSH= (M**2*PI*FKC*L0K) / ( SIGH AK*DDS) 

C 

C**** A RHA 1U RE LCSSES 

PA= (AJA**2/ (SIGHAA*SFA)) *THBAE*RAO**2* (1.0-X**2) *LOA*NPA 
C 

C**** FIELD AT SHIELD RADIUS 
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BRS=MU* ((4.0* A JP*SC 11) / ( BB*PI) ) *HS*Z** EB* (1.0-T**BB) 

C R SO 15 SHI 11C OUTER RADIOS 

C SHIS ID IS DESIGNED FOR UKIEORM FLUX CENSITT 

RSO* RS* (1. C*PR3/ (BSH A X*P) ) 

C 

€•*** PI ELD AT AN INNER CORNER OP FIELD WINDING 
HR =0 
HTH *C 

C ESTIMATE PIELD INTENSITY AT INNER RADIOS AND HIGHEST ANGOLAB 

C EXTENT OP THE PIELD WINDING, HARMONICS 1 TO 19 AND BOTH RADIAL AND 

C AZIMUTHAL COUPON ENTS ARE INCLUDED. 

DO 1500 1=1,20,2 
F=PLG AT ( I) *P 
G=FLCAT (I) *THW EE/2. 0 

1530 H (I) =2 . 0* A JE* SIN (G)/ (FLOAT (I)*PI* (2.0-F) ) *R FI*Y** (F-2. 0) 

1 *CM (P,Y,Z) * (4 . 0-P* * 2) *E 
DO 16C0K=1 ,20, 2 
EE=P10 AT ( PC ) ’THWPE/2.0 
fl R=H B*H (K) * SI N (EE) 

1600 HTH=HTH*H (K) *CCS (EE) 

HMAX-SQRT(HR**2*HTH**2) 

C 

C***‘ ROTOR CRITICAL SPEED 

C COMPUTE ROTOR CRITICAL SPEED USING A SIMPLE BENDING MOMENT MODEL 

C TORQUE TUBE STIFFNESS ONLY IS USED 

C MASS HA INCLUDES TORQUE TUBE, SHIELD, AND PIELD 

M A= (ROSS* (RFI**2-RFI*-2) *R CCU* SEP* RF 0*2* (1 . 0-Y**2) 

1 *R OAL* (FKQ**2-FKI**2) ) *PI*FOE* ( (R PO* GF K- GKI) **2-RFO**2) *PI 
C IN IS TORQUE TUBE FOMENT CP INERTIA 

IN - (3 • 14 16/4. 0) * (FEI**4-RBI**4) 

OH GCR T= 9. 875*SQF T ( E*IN/ (LPF**4*HA) ) 

OHBPM=CMGC RT/ ( 2* 0* FI) *60.0 
C 

C**** WEIGHTS OP FAJCR MATERIALS 

C STAINLESS STEEL (TORQUE TUBE), COPPER (ARMATURE) , STEEL (IRON SHIELD) 

C AND AIUEI SUM (DAMFEP) 

M SS- PI * (R-I**2-RPI**2) *LBP*ROSS 
HK=PI*R KO**2* ( 1. 0-ZZ**2) *LCK*FCAL 
HA=PI*FAO-*2* ( 1. 0-X**2) *LO A*ROCU*SP A 
MS=PI* (RSO **2“PS**2) *LCA-FCFE 
C MB IS WEIGHT OE BINDING MATERIAL 

HR=PI* ( (RPO^GF K-GKI) **2-RPC**2) *ICP 
C 

C**** SU PERCON DUCTIO R REQUIREMENTS POB COST ESTIMATES 
C RESULT IS IN A FPER E-TU R NS-FETERS 

ATM* (AJF*THWFE-R PO**2* ( 1 . 0-Y** 2) /2. 0) * 2. 0*LOF 
C 

C**** STATCR CORE LOSSES 

C PPKG IS CORE LOSS DENSITY IN WATTS/KG 

PPKG=PZ* (BRS/BSHAX)**GABMA 
PCOR E* M 5*PFKG 
C 

C** # * COST FUKCTICN 

C C IS COST OF MACHINE IN WEIGHT 

C*HSS*FA*HK*ES*HB 
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ALOS=P A *PCOR F*PS H 
C V LOS IS IN KG FEB WATT 
VLOS=C/ALOS 
W ATLC5 =ALOS 
KGLOS*tf LOS*ALOS 
PFLS=1.0 
C 

C**** PEN ALTY FUNCTION FCR FI ILL CURRENT LI BIT 

C FIELD CHARACTERISTICS A FE IHPUI AS FIVE POINTS IN H-J PLANES 

C H IS IN HC ARRAY, J IS IN JC ARRAY 

C HIS ASSUMED CONSTANT IN A SWING 

C FOR A FIXED H, LINEAR INTERPOLATION IS USED TO OBTAIN CRITICAL 

C CURRENT DENSITY 

C HI AX IS MAXIMUM FIELD INTENSITY 

IF (UMAX .LE. HC ( 1) ) GO TO 10 
PFFC= 1 0 **4 0 
GO TC 11 

10 DO 1700 1=2,5 
K*I-1 

L=I ♦ 1 

IF (Hi AX .GT. HC ( I) .AND. HM AX . LE. HC ( K) ) 

1 JSR= JC (I) ♦ (JC(l) - JC (I) ) * (HPAX-HC(I) ) /(HC(L) -HC(I) ) 

1700 CONTINUE 

IF(RHAX .LT. HC ( 5) ) JSR = JC(5) 

PFPC=.9*.1*(AJF*IF1/JSF) ** 15 

c 

c »«*« PENALTY FUNCTION FCR SHAFT STRESS 

11 PFSI = . 9*. 1* (TD/TMAX) **15 
C 

C**«* PENALTY FUNCTION FOR SHAFT CRITICAL SPEED 
C PROVISION IS HADE TO FCFCE LCW CRITICA SPEED AND TO KEEP 

C CRITICAL SFEFD AWAY FROM OPERATIN SPEED 

PFCS= . 9*. 1* (OH EG A/ (P-OMGCRT) ) **3* . C01* ( P* OMGCRT/ ( P* OMGCRT- 
1 0flFG A) ) **2 
C 

PENALTY FUNCTION FCR SHIELD FLUX LIMIT 
PFSF= , 9 ♦ . 1 * (BPS/BSKAX) **15 
C 

C**** PENALTY FUNCTION FOR DAMPER STRESS 
PFDS=.9*. 1* (STB/DHAX) ** 5 
C 

C**** PENALTY FT7NCTICN FCR ARMATURE INSOLATION AND DIAMETER 
PFAI = .9*. 1*(GAI/(BAI-RKO) ) **15 
PF AC= . 9+ . 1* (GAI/ (RS-RAC) ) **15 
C 

C**** FINAL COST FUNCTION 

CF=C-PFSI*PFCS*PFSF*PFFC*PFDS*PFLS*PFAI*PFAC 

IF (NUM .GT. 0) GO TO 100 

V 1 1= V A/1E06 

V12=HMAX*MU 

WRITE (6,701) 

WRITE (6,703) 

WRITE (6,7T5) 

WRITE (6,707) HP 
WRITE (6,709) VII 
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29 1 

292 

29 3 

294 

295 

296 

297 
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299 

300 
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30P 

309 
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312 
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317 
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WRITE (6, 


r 7 1 1 ) 


PF 


WRITE (6, 


,713) 


RPI1 


WRITE {6, 


,715) 


P 


WRITE {6 


,717) 


VI 


WRIir £6, 


,716) 


A JA 


WRITE (6 , 


,719) 


KVAFH 


WRITE (6, 


,70 3) 




WRITE (6, 


,721) 




WRIIE (6 , 


,723) 


V £6) 


WRITE (6 


,725) 


v £7) 


WRITE {6, 


,727) 


RAO 


WRITE (6, 


, 729) 


LCA 


WRITE {6 ( 


,731) 


LA 


WRITE (6, 


,73 3) 


LAS 


WRITE £6 ( 


,735) 


V FT 


WRITE £6, 


,737) 


N A AT 


WRITE £6 


,739) 


SFA 


WRITE £6 


,741) 


NPA 


WRITE £6 


,743) 


SIGNAA 


WRITE £6 


,74 5) 


THWAE 


WRITE (6 


,703) 




WRITE £6 


,751) 




WRITE £6 


,753) 


V £2) 


WR ITE £6 


,755) 


V £1) 


WRITE £6 


,757) 


V £3) 


WRITE £6 


, 759) 


LA 


WRITE £6 


,76 1) 


LCF 


WRITE £6 


,763) 


A JR 


WRITE £6 


,765) 


SFP 


WRITE £6 


,767) 


THW*E 


WR ITE £6 


,769) 


NEAT 


WRIT? |6 


,771) 


I F 1 


WR ITE £6 


,773) 


V 12 


WRITE (6 


,775) 


XD 


WRITE £6 


,777) 


XEP 


WRITF £6 


,779) 


XCPP 


WRITE £6 


,703) 




WPITE £6 


,701) 




WRITE £6 


,783) 


V£4) 


WRITE (6 


,78 5) 


RKO 


WR ITE £6 


, 767) 


v £5) 


WRITE (6 


,709) 


LCK 


WRITE (6 


,791) 


LAK 


WRITS (6 


,79 3) 


SIGHAK 


WRITE £6 


, 703) 




WPITF £6 


,706) 




WRITE £6 


,797) 


XT 


WRITE £6 


,796) 


XI 


WRITE £6 


,799) 


X 2 


WRITE £6 


,703) 




WRITE £6 


,001) 




WRITE £6 


,003) 


OflGCRT 


WRITE (6 


, 005) 


LPR 


W RITE £6 


,703) 




WRITE £6 


,811) 
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345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

36 0 

361 

36 2 

363 

36 4 

365 

366 

367 

366 

369 

370 

371 

372 

37 3 

374 

375 

376 

377 

378 

379 

380 

301 

382 

38 3 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

39 4 

395 

396 

397 



WRITE 


(6 #8 1 3) 


MSS 


WRITE 


(6,815) 


HK 


WRITE 


{6,817) 


MA 


WRITE 


{6,819) 


fie 


WPITE 


(6,821) 


MS 


WRITE 


(6,703) 




WF ITE 


(6,831) 




WRITE 


{6,03 3) 


BFS 


WRITE 


(6,035) 


RSO 


WRITE 


(6,837) 


RS 


WHITE 


(6,83 9) 


BSK AX 


WHITE 


(6,70 3) 




WRITE 


(6,841) 




WRITE 


(6,843) 


PA 


WRITE 


(6,945) 


PCORE 


WRITE 


(6,84 7) 


PSH 


WRITE 


(6,849) 


WATIOS 


WRITE 


(6,851) 


WLOS 


WRITE 


(6,853) 


KGLCS 


WRITE 


(6,703) 




WRITE 


(6,671) 




WPITF 


(6 ,873) 


FCCO 


WRITF 


(6,675) 


FCFI 


WRITE 


(6,877) 


RCAL 


WRITE 


(6,879) 


FOSS 


WR ITE 


(6, BBC) 


RCB 


WRITE 


(6 ,703) 




WHITE 


(6,881) 




WPIT2 


(6,983) 


PFSI 


WRITE 


(6,885) 


PFCS 


WRITE 


(6,887) 


PFSF 


WRITE 


(6,889) 


PFFC 


WRITE 


(6,991) 


PFCS 


WRITE 


(6,89 0) 


PFAI 


WRITE 


(6,892) 


PFAC 


WRITE 


(6,703) 




WRITE 


(6,901) 




WRITE 


(6,905) 


C 


WPITE 


(6,907) 


CF 


WRITE 


(6,703) 




W RIT? 


(6,911) 




WRITE 


(6,913) 


DM A X 


WRITE 


(6,915) 


THAI 


WRITE 


(6,917) 


E 


WRITE 


(6,919) 


PR 


WRITE 


(6 ,703) 




WRITE 


(6,5000) 



701 FORMAT {131 ,10X ,40H SOPERCOHDUCTING G EH ES ATOR /KOTOR DESIGH 
703 FORMAT (/72H 



1 



705 


FOPMAT 


(13H 


* * 


707 


FORMAT 


(5X, 


46 H 


7C9 


FORHAT 


(5X. 


U6H 


711 


FOR K AT 


(SX, 


468 


713 


FORKAT 


(SX. 


46H 



— "/) 

RATING **) 

FATED POWER (8 F) 

FATED POWER (KVA) ..... 

POWER FACTOR 

MECHANICAL SPEED (BPK) 



F10.0) 
>F1 C.O) 
F 10. 3) 
,F10.0) 
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398 

399 

400 

401 

402 

403 

404 

405 

406 

4 07 

40 8 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

4 20 

421 

422 

423 

424 

425 

4 26 

427 

428 

429 

430 

431 

432 

433 

434 

435 

4 36 

437 

4 38 

439 

44C 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 



715 FORMAT (5X, 46H NUMBER OF POLE PAIBS # F10.0 

717 FORMAT (5X # 46R TERMINAL VOLTAGE (V) ,F10.0 

718 POR MAT ( 5X* 46 R ARMATURE CURRENT (A) E10.4 

719 FORMAT (5X,46H PER UNIT POWER RATING (P.U.J ,F10.2 

721 FORM AT ( 1 5H •• ARHATUPE **J 

723 FOR MAT (5X , 46H ARMATURE THICKNESS (MJ ,F1C.4 

725 FORMAT (5X f 46H ARMATURE TC CCRE GAP { MJ ,F10.4 

727 PO RMAT (5X f 46H ARMATUPE OUTER RADIUS (MJ ,F10.4 

729 PORHAT (5X,46H AVER ALL ARMATURE LENGTH (MJ ,F10.4 

731 FORMAT (5Y, 46H ACTIVF ARMATURE LENGTH (H) ,F10.4 

733 PORHAT (5X,46H STRAIGHT SECTION LENGTH (MJ ,P10.4 

735 FORMAT (SX, 4611 VOLT PFR TURN (RMS) ,P10.5 

737 FORMAT (5X,46H ARMATURE AMPERE-TURNS (FMS/PHASE) ,E10.4 

739 FORMAT (5X # 46H ARMATURE WINDING SPACE PACTCF ,F10.4 

741 FORMAT (5X,46H NO. OF ARMATURE PHASES ,F10.4 

743 FORMAT (5X,46H ARMATURE CONDUCTIVITY (S/B) ,E10.4 

745 FORMAT (5X, 46H ARMATURE AAGLE (RAD) P10.4 

751 FO RMAT (23R FIELD WINDING **J 

753 FORMAT (5X, 46H FI EL E THICKNESS (M) # F10.4 

755 FORMAT (5X # 46H FIELD INNER RADIUS (H) # F10.4 

757 FORM AT (SX # 46H FIELD TO DAMPER GAP (ft) # F10.4 

759 FORMAT (5X,U6H ACTIVE MACHINE LENGTH (MJ # F10.4 

761 FORMAT (5X, 46H OVERALL FIELD LENGTH (MJ ,F10.4 

763 FORMAT (5X,46l! FIELD CURRENT DENSITY (A/M**2) # E10.4 

765 FORM AT (5X,46U FIELD WINDING SPACE FACTOR # F10.4 

767 FORMAT (5X,46H FIELt ELECTRICAL WINCING ANGLE (RAD) # P10.4 

769 FORMAT (5X f 46H FIELD AMPERE-TURNS (A-T) ,E1C.4 

771 FORM AT (5X, 46H MAX PER UNIT PIELC CORBEKT (P.U.J ,P10.2 

773 FORMAT (5X.46H MAXIMUM FIELD (TESLA) # F10.0 

775 FORMAT (5X, U6H SYNCHRONOUS REACTANCE # F10.4 

777 PORMAT (5X,46H TRANSIENT REACTANCE ,F10.4 

779 FORM AT (5 X # 46 H SUfiTR A NS I E NT REACTANCE ,F10.4 

781 FORMAT (13M ** DAMPER **) 

783 FORM AT (5 T , 46H CAMP FR THICKNESS (H) # F10.4 

785 FORMAT (5X,46H DAMPER OUTER RADIUS (MJ # P10.4 

787 FORMAT (5X,46H CAMPER TO ARMATURE GAP (MJ ,F10.4 

789 FORMAT (5X,46H OVERALL CAMPER LENGTH (MJ # F10.4 

791 FORMAT (5X,46H ARMATURE COUPLING LENGTH (M) ,P10.4 

793 FORM AT (5X # 46H CAHPFR CONDUCTIVITY ,E10.4 

796 FORMAT (16U ** STABILITY **) 

797 POR MAT (5X,46H TRANSFORMER REACTANCE f F10.4 

798 FORMAT (5X,46H REACTANCE UNFAULTED LINE # F10.4 

7Q9 FORMAT (5X,46H READTANCE PAULTED LINE ,F10.4 

801 FORMAT (24H ** NATURAL FREQUENCY **) 

803 FORMAT (5X # 46,H ROTOP CRITICAL SPEED (RPM) # F10.1 

805 FORMAT(5X,46H BEARING SPAN (MJ ,F10.4 

811 FORM *T ( 1 3H *+ WEIGHT **J 

813 FCPMAT (5X,U6H STAINLESS STEEL SUPPORT (KG) # F10.4 

815 FCR;1AT(5X,U6H SHIELD WINDING (KG) P10.4 

817 FORMAT (5X,46H ARMATURE (KG) F10.4 

819 FORMAT (5X.46H PIN DING MATERIAL (KG) ... # P10.4 

821 FORMAT (5X, 46H STATCR CCRE (KG) # F10.4 

831 FORMAT (26H ** FEROHAGNETIC SHIELD **J 

833 FORMAT (5X,46H FLUX AT SHIELD RADIUS (TESLA) # F10.4 

835 FORMAT (5X # 46 H SHIELD OUTER RADIUS (MJ ,F10.4 
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453 

454 

455 

456 

457 

450 

454 

460 

46 1 

462 

463 

464 

465 

466 

467 

460 

469 

470 

471 

472 

473 

474 

475 

476 

477 

470 

479 

400 

401 

402 

403 

404 

405 

406 



037 romT(5X # 46H SHIELD INNER RADIOS (R) ,M0.4) 

039 POHM *T(5X,46H PAX. SIIIEID PLUX DENSITY (TESLA) ,P10.4) 

041 FOR RAT ( 13!l ♦* LOSSES «♦) 

043 PORRAT (5X,46H ARM A TOR E LOSSES (WATTS) # E10.4) 

045 P03HAT (5X,46H STATOR CORE LOSSFS (WATTS) # P10.4) 

047 PORT AT (5X # 46H NEGATIVE SEQUFWCF LOSSES (WATTS) # E10.4) 

049 PORRAT (5X,46H TOTAL LOSSES (WATTS) # E10.4) 

051 FORMAT (5X,46H COST OF LOSSES (KG/WATT LOST) ,F10.4) 

053 FORMAT (5X.46H TOT A l COST CP LCSSES (KG) ,E10.4) 

071 FORMAT (16H ** TENSITIES **) 

073 PORRAT (5X,46fl COPPFP (KG/K**3) ,F10.4) 

075 PORRAT (5X, UGH I PON (KG/R**3) # F10.4) 

877 FORMAT (5X,46R ALUHINUR (KG/r>*3) # P10.4)> 

079 F ORRAT (5X , 46H STAINLESS STEEL (KG/fl**3) # F10.4) 

000 PORI AT (5X.46H PIN DING RATFFIAL (KG/R**3) ,P10.4) 

001 FORMAT (74H ** PENALTY FUNCTIONS **) 

003 PORT AT (5X, 46H SHAFT STRESS ,P10,4) 

085 FORMAT (5X.46H SHAFT CRITICAL SPEED ,P10.4) 

087 FORRAT (5X,46H SHIEID FIUX LIHIT # F10.4) 

0B9 Foa RAT (5X, 46H FIELD CURRENT LIMIT # P10.4) 

091 FORRAT (5X.46H EAKPEP STRESS ,P10.4) 

090 FOR RAT (5X,46H ARMATURE INSULATION THICKNESS # F10.4) 

092 FORRAT (5X # 46H ARMATURE DIAMETER P 1 C . 4 ) 

901 FORRAT (11H ** COST ♦*) 

905 FORRAT (5X # 46H COST FUNCTION ,E10.4) 

907 FORM AT (5X, UGH PENA IIZEE CCST FUNCTION # E10.4) 

911 FORM AT ( 25H ** MATERIAL CONSTANTS *•) 

913 FO p M AT (5X # 46H MAX SHEAR STRESS IN DA HPER MATERIAL # E10.4) 

915 FORMAT (5X,46H RAX SHEA B STRESS IN TOFQIE TUBE MATERIAL .... # E10.4) 

917 FO RR AT (5X # 46H YOUNG* S MODULUS E 1 0. 4) 

919 FORMAT (5X # 46H POISSON*S RATIO # P10.4) 

5000 FO RM AT ( * 1 • ) 

100 RETURN 
FNC 
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487 




FUNCTION CM(P,X,W) 




c**** 


CALC GEOMETRIC COEF. CM **** 


488 




IF ( P-2 . 0 ) 1000 , 1100 , 1000 


489 


1100 


CM=0 . 125* ( -AL0G10(X )+0. 03125* (1.0-X*X)*W**4) 


490 




GO TC 1111 


491 


1000 


AA=2.C-P 


492 




BB=2 . 0+P 


493 




CC=2.0*P 


494 




CM= ( ( 1 . 0-X** AA)+ ( AA/BB)* ( 1 . 0-X** BB)* (W**CC ) ) 
l/(P*(4.0-P**2)) 


495 


mi 


RETURN 


496 




END 


497 




FUNCTION CS(P,X,W) 




c**** 


CALC GEOMETRIC COEF. CS **** 


498 




IF (-2.0) 1000,1100,1000 

CS= ( ( X** 4* ( AL0G10 (X) ) ) /20 . 0 ) + ( ( 1 . 0-X** 4/8 . 0) 
l+( ( ( (1.0-X**4)**2)/l6.0)*W**4) 


499 


1100 


500 




GO TO 1111 


501 


1000 


AA=2 . 0-P 


502 




BB=2.0+P 


503 




CC=2.0*P 


504 




C S= ( ( AA- ( 4 . 0* X** BB)+ ( BB* X** 4 ) + ( 2 . 0* AA/BB* 

1 ( ( 1 . 0-X** BB)** 2) )* (W**CC ) )/(P* (4.0-P**2) ) ) 


505 


ini 


RETURN 


50 6 




END 
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APPENDIX B 



Output of 20,000 Horsepower Generator Optimization 
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superconducting generator/mctor design 



*• RATING ** 

RATED PGwER (HP) 2DK7. 

RATED PGWER (MV A) 15. 

PCWER FACTOR l.JCO 

MECHANICAL SPEED (RPM) 169D. 

NUMBER DF POLE PAIRS 1. 

TERMINAL VOLTAGF IV) 1. 

ARMATURE CURRENT IA) 0.35C0E 07 

PER UNIT POWER RATING (P.U.) . 1.00 



*♦ ARMATURE ** 

ARMATURE THICKVESS I M ) 0.0957 

ARMATURE TC CORE CAP IN) 0.0272 

ARMATURE CUTER RACIUS (M) 0.3058 

AVER ALL ARMATURE LENGTH (H) 1.0R7B 

ACTIVE ARMATURE LENGTH (M) 0.6846 

STRAIGHT SECTION LENGTH I M ) 0.5740 

VOLT PER TURN (RMS) 45.4^466 

ARMATURE AMPERE-TURNS (RMS/PHASE) 0.918PE C5 

ARMATURE WINDING SPACE FACTOR 0.3000 

NC. 3F ARMATURE PHASES 3.OJ0O 

ARMATURE CONDUCTIVITY (S/M) D. 60002 Vb 

ARMATURE ANGLE (RAD) 1.04 7D 



FIELD WINDING ** 

FIELD THICKNESS I M ) C.026O 

FIELD INNER RADIUS (M) 0.0°72 

FICLD TO DAMPER GAP I M ) 0.0225 

ACTIVE MACHINE LENGTH (M) 0.6846 

CVCRALL FIELC LENGTH I M ) D.9059 

FIELD CURRENT DENSITY (A/K**2) 0.1575E C9 

FIELD WINCING SPACE FACTOR D.50C0 

FIELD ELECTRICAL WINDING ANGLF (RAD) 2. 0940 

FIELD AMPERE-TURNS I A - T ) D.9819E 06 

MAX PER UNIT FIELC CURRENT (P.U.) 1 . 10 

MAXIMUM FIFLC (TFSLA) 4. 

SYNCHRONOUS REACTANCE 1.38U0 

TRANSIENT REACTANCE 1.1845 

SUBTRANSIENT REACTANCE 1.2017 



*• CAMPER ** 

DAMPER THICKNESS (M) 0.^48 

DAMPER OUTER RADIUS (M) . 0.1914 

CAMPER TO ARMATURC GAP (M) 0.0227 

Overall dampf.r length <mj 1.0226 

ARMATURE COUPLING LENGTH IM 1 0.<U19 

DAMPER CONDUCTIVITY 0.2000F OP 



•* STABILITY •* 

TRANSFORMER REACTANCE 0.1000 

REACTANCE UNTAULTED LINE 0.10C0 

READTANCE FAULTED LINE 0.15C0 



• * NATURAL FREQUENCY ** 

RCTOR CRITICAL SPEED (RPM) 5P6.3 

SEARING SPAN (M> 1.5226 



** WEIGHT ** 

STAINLESS STEEL SUPPORT (KG) 68.69C0 

SHIELD WINDING (KG) 126.4667 

ARMATURE (KG) 456.5620 

BINDING MATERIAL (KG) 0.001/ 

STATOR CORE (KG) 2523.6330 



•* FERCMAGNET IC SHIELD ** 

FLUX AT SHIELD RACIUS (TESLA) 0.6361 

SHIELD OUTER RADIUS ( M ) 0.4595 

SHIELD INNER RADIUS ( M ) 0.3)70 

MAX. SHIELD FLUX CENSITY (TESLA) 1.7500 



•• LOSSES ** 

ARMATURE LOSSES (WATTS) 0.U77E C6 

STATOR CORE LOSSES (WATTS) 599.4)70 
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NEGATIVE SEQUENCE LOSSES < WATTSI 0.3569F 03 

TCTAL LOSSES I WATTS I **♦* 0.U86F 06 

CCST OF LOSSES (KC/WATT LCST ) 0.0268 

tCTAL COST OF LOSSES (KG) 0.3175F 0* 



** OENSITIES *♦ 

COPPER (KG/M**3 ) 8800.0000 

IRON ( KG/M**3 ) 7500.3000 

ALUMINUM (KO/M**3) 2600.0000 

STAINLESS STEEL <KG/M**3) fiOOO.OOOO 

BINDING MATERIAL (KG/M**1) 1800. OOCO 



PENALTY FUNCTIONS ** 

ShAFT STRESS 1.0330 

SHAFT CRITICAL SPEED 0.9344 

SHIELD FLUX LIMIT 0.0000 

FIELD CURRENT LIMIT 0.9051 

OAMPER STRESS 0.90CG 

ARMATURE INSULATION THICKNESS 0.9152 

ARMATURE DIAMETER 0.9010 



*«■ CCST ** 

CCST FUNCTION 0.3175E 04 

PENALIZED COST FUNCTION ••••• ....0.1794E 04 



MATERIAL CONSTANTS +* 

MAX SHEAR STRESS IN DAMPER MATERIAL 0.2400E 39 

MAX SHEAR STRESS IN TORQUE TUOE MATERIAL ....0.4530E C9 

YOUNG'S MODULUS 0.2000C 12 

PCISSON'S RATIO 0.30^0 
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APPENDIX C 



Output of 40,000 Horsepower Motor Optimization 



SUPERCONDUCTING GENERA TOR /MOTOR OESIGN 



•• RATING ** 

RATED POWER I Hf> ) 40214. 

RATEO POWER IMvA) 30. 

POWER FACTOR l.OCO 

MECHANICAL SPEED IRPM) 20C. 

number OF POLE PAIRS 3. 

TERMINAL VOLTAGE <V) 1. 

armature current i a ) o.3sooe or 

PER UNIT POWER RATING IP.U.) 1.00 



•• ARMATURE ** 

ARMATURE THICKNESS IM 1 0.1739 

armature tc coir gap imj 0.0309 

ARMATURE OUTER RADIUS IM) 0.7292 

AVER ALL ARMATURE LENGTH |M) 1.7H77 

ACTIVE ARMATURE LENGTH |M) 1.5085 

STRAIGHT S r CT I ON LENGTH I M > 1.3596 

VCLT PER TURN IRMS)... * 19.49800 

ARMATURE AMPERE-TURNS IRMS/PHASE) 0.4092E C6 

ARMATURE WINDING SPACE FACTOR 0.30C0 

NC. OF ARMATURE PHASES 3.0000 

ARMATURE CONDUCTIVITY I S/M) 0.6000E 08 

ARMATURE ANGLE IRAQ) 1.0470 



•• FICLO WINDING ** 

FIELO THICKNESS I M) 0.0310 

FIELO INNER RADIUS IM) 0.4315 

FIELD TO DAMPER GAP IM) 0.0225 

ACTIVE MACHINE LENGTH IM) 1.5085 

OVERALL FIELD LENGTH IM) 1.8065 

FIELO CURRENT DENSITY IA/M**2) 0.1432E 09 

FIELO WINCING SPACE FACTOR 0.5000 

FIELO ELECTRICAL WINDING ANGLE IRAD) 2. 3940 

FIELD AMPERE-TURNS I A-T I 0.4152E 07 

MAX PER UNIT FIELD CURRENT IP.U.) 1.14 

MAXIMUM FIELD (TESLA) 4. 

SYNCHRONOUS REACTANCE 1.0C49 

.TRANSIENT REACTANCE 0.8668 

SUBTRANSIENT REACTANCE 0.8643 
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*• DAPPER •• 

CAPPER THICKNESS (M ) 0.0449 

DAMPER OUTER RADIUS (H) 0.5298 

CAMPER TO ARMATURE CAP CM) 0.0255 

OVERALL DAMPER LENGTH IK) 1.8468 

ARMATURE CCU°LING LENGTH (M) 1.6978 

CAMPER CONDUCTIVITY 0.20C0E 08 



•A STABILITY ** 

TRANSFORMER REACTANCE 0.10C0 

REACTANCE UNFAULTEO LINE 0.13C0 

REACTANCE FAULTED LINE 0.1000 



** NATURAL FREQUENCY ** 

RCTOR CRITICAL SPEED I RPM ) 1905.9 

8EARING SPAN ( M ) 2.3468 



•* WEIGHT ** 

STAINLESS STEEL SUPPORT ( KG) 2036. 2830 

SHIELD WINDING (KG) 687.1162 

ARMATURE (KG) 3311.5350 

BINDING MATERIAL (KG) 0.0132 

STATOR CORE (KG) 7309.2100 



• • FERC’TAGNETIC SHIELD ** 

FLUX AT SHIELD RADIUS (TESLA) .., 

SHIELD OUTER RADIUS ( M ) 

SHIELD INNER RADIUS ( M ) 

MAX* SHIELD FLUX DENSITY (TESLA) 

LOSSES 

ARMATURE LOSSES (WATTS) 

STATOR CORE LOSSES (WATTS) 



0.7367 
D . 8667 
0.76CI 
1.7500 



D. 8535C P 6 
2428. 802C 
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